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adapt  conputer- 


technique.  Tha  analytical  approach  is  tc 
aidsd  dasign  cptimiaaclon  to  find  optinuD  paraDstara 
automatically. 

Thraa-diaenaional  approxiaata  dose  models  ara  davaloped 
to  aimulate  the  exact  dose  aodal  using  a spherical  or 
cylindrical  coordinate  system.  Optimum  parameters  are  found 
much  faster  with  the  use  of  computer-aided  design 

The  iaplementation  of  computer-aidad  dasign  algorithms 
with  the  approximate  dose  model  and  the  application  of  the 
algorithms  to  several  cases  are  discussed.  It  is  shown 
that  the  approximate  dose  model  gives  dose  distributions 
similar  to  those  of  the  exact  dose  model,  which  maXas  tha 
approximate  dose  model  an  attractive  alternative  to  the 
exact  dose  model,  and  much  more  efficient  in  terns  of 
computar-aided  design  and  visual  optimisation. 


Staraotactic  radioaurgary  is  a technics  for 
obliterating  intracranial  tuaors  which  are  inaccessible  or 
unsuitable  for  open  surgical  technigues,  using  tightly 
colligated  beams  of  ionizing  radiation.  The  aim  of 
stereotactic  radiosurgery  la  to  deliver,  with  a high  degree 
of  spatial  accuracy,  a large  radiation  dose  to  the  target 
volume  within  the  brain,  while  maintaining  the  smallest 
passible  dose  to  the  remainder  of  the  brain  tissue.  At 
present,  there  are  three  types  of  radiation  beans  that  can 


be  used  for  radiosurgery:  focused  beans  obtained  from 
specially  designed  units  incorporating  several  hundred 
cobalt  sources  (the  eamma  Knife)  (Larsson  et  al.,  1974: 
□ahlin  & Sarby,  1975:  LaKsell,  1983):  heavy  charged  particle 
beams  (protons,  deuterons,  helium  ions,  etc.)  obtained  from 
cyclotrons  or  synchrocyclotrons  (Lyman  8 Howard,  1977:  Lyman 
et  al.,  1986):  and  x-rays  obtained  from  medium  energy  (4  to 
10  MV)  isocentrlcally  mounted  linear  accelerators  (tlMAC) 
(Helftz  et  al.,  1984:  Colombo  at  al.,  19BS;  Hartmann  et  al.. 


1988:  Friedman  S Bova,  1989).  LtNAC-based  stereotactic 


radiosurgery 


expensive  alternative 


proven  radioaurgical  technlgues.  The 


on  the  conbinaClon  of  multiple  isocentric  arc  irradiation 
with  small  fields  centered  in  the  stereotactic  target.  This 
optimizes  the  dose  falloff  outside  the  target  volume. 

Radiosurglcal  treatment  of  the  brain  requires  a 
detailed  thrae-dimensicnal  (3-D)  treatment  planning  system, 
and  intervening  and  surrounding  tissues  must  be  protected 


i irradiation  while  a high  d 
to  the  shape  of  the  target  is  produced, 
optimal  radiosurgery  planning  system  which  uses  3-D  patient 
data  and  treatment  parameters  repreeante  a significant 
challenge.  This  is  in  part  due  to  the  lacJc  of  3-D 
information  about  target  volumes  and  anatomic  structures, 
and  also  by  the  many  beam  parameters  involved  i 


planning. 


with  the  availability  of  computed  tomography  fCT]  and 
magnetic  resonance  (MR)  Imaging  devices,  the  ability  to 
visualize  anatomic  structures  has  been  greatly  enhanced. 
Many  new  tools  have  been  developed  end  evaluated  which 
permit  the  generation  of  3-D  displays  of  these  structures, 
the  accurate  specification  of  bean  direction  and  size,  the 
display  of  dose  distribution  in  three  dimensions,  and  the 
eveluation  at  treatments  for  optimization.  These 
improvements,  however,  have  also  added  to  the  complexity  of 
evaluating  a given  plan.  The  full  3-D  extent  of  the  tumor, 
the  neighboring  anatomy,  their  relationship  to  the  planned 


radiation  fields,  and  the  resulting  dose  distribution  must 
all  be  considered. 

The  use  of  graphic  displays  to  present  the  information 
obtained  from  ioproved  3-D  planning  can  be  helpful  in 
acquiring  an  overall  impression  of  tha  adequacy  of  a 
particular  plan.  However,  3-D  planning  uses  a large  volume 
of  patient  data  and  many  beam  parameters,  which  nay  vary 
with  the  treatment.  Furthermore,  since  3-D  calculations 
tahe  a very  long  time,  interactive  modification  of  treatment 
using  graphic  displays  nay  not  be  practical.  Another 
alternative  to  complex  3-D  optimised  treatment  is  to  use  an 
analytical  method  of  handling  treatment  plans 
quantitatively.  Sufficient  information  on  biologic  dose 
response,  tumor  control  probability,  and  normal  structure 
complication  rate  is  needed  to  quantify  optimal  treatments. 
Since  no  detailed  information  about  biologic  optimization  is 
available  at  present,  the  evaluation  method  is  based  on 
physically  optimizing  dose  distribution. 

A possible  solution  for  3-D  treatment  plan  optimization 
is  to  utilize  computer-aided  design  optimization  techniques 
with  proper  objective  functions  to  represent  the  physical 
optimization  criteria.  The  use  of  analytic  formalism  as  an 
objective  function  for  automatic  optimization  has  been  tried 
with  little  clinical  success  and  is  still  being  studied. 

This  effort,  however,  was  limited  to  optimizing  simple 
linear  variables  such  as  beam  weights  or  treatment  times, 


reflected  only  two-dinensional  consldecatlone. 


based  radiosurgery  uses  many  noncoplanar 
evaluation  technigue.  Accordingly  aany 
bean  paraneters  and  complex  3-D  calculation  procedures 
:he  dose  optinization.  Although  the  most 


difficult  problem  ie  t 

radiosurgery  may  not  b 
visual  and  computer-ai 
The  aim  of  this  w 
of  planning  optimal  do 


at  of  specifying  the  optimization 
of  3-D  treatment  planning  for 
possible  without  the  use  of  both 
ed  design  optimization, 
rk  is  to  investigate  the  feasibility 
e distributions  for  LINAC-based 


stereotactic  radiosurgery. 


both  experimental  optimization  with  visual  evaluation  of 
dose  distribution,  and  analytical  optimization  through 
computer-aided  design  optimization.  The  experimental 
approach  is  based  on  ehaping  the  target  volumes  through  the 
use  of  multiple  isocenters  or  field  shaping  techniques.  The 
analytical  approach  is  to  adapt  computer-aided  design 

3-D  dose  algorithm  to  avoid  dose  to  critical  organs  and  find 
optimum  arc  positions  automatically  while  maintaining  a dose 
distribution  shaped  to  the  target. 


investigated  some 
3-0  optimization. 


associated  with  e comprehensive 


Patient  Data 


The  final  set  of  patient  data  is  represented  by  a 
series  of  transverse  sectional  contours  along  the  patient. 
This  set  of  3-D  patient  data  contains  external  contours, 
target  cross-sections,  and  relevant  normal  tissue  outlines. 

One  difficulty  in  obtaining  an  optimal  dose 
distribution  is  due  to  different  patient  head  geonetrlaa  and 
target  positions.  Another  problem  originates  from  the 
difficulty  of  outlining  the  3-D  solid  shape  of  Che  target 


and  critical  organs.  It  is  inconvenient  to  find  the  optimum 
parameters  for  each  patient  condition  every  time.  A better 
approach  is  to  use  a general  simplified  patient  model  which 
can  simulate  accurately  the  real  head  shape  or  target  shape 
without  any  significant  difference  between  the  simplified 
model  and  real  patient  geometry. 

A spherical  head  model  with  a centered  target  poeiticn 


is  assumed  for  the  experimental  optimisation  study,  since  it 
is  geometrically  simple,  and  represents  a useful  result  with 
a RAHDO  head  model.  The  dose  shape  is  not  much  changed 
except  in  an  extrema  casa  fnear  the  surface)  as  the  target 
position  varies.  Therefore,  useful  information  such  as  the 
isodose  shape,  or  target  margin,  obtained  from  this 


reference  model  can  be  used  as  a guideline  for  an  individual 


patient  model. 


The  next  prohles  is  to  identify  a 3-D  target  shape  and 
the  critical  organs.  For  sinplicity,  the  target  shape  and 
critical  organs  are  approxiaated  by  ideal  geonetries  (e.g. 
cylindrical,  spherical,  cone,  etc.),  which  sufficiently 
cover  and  siaulate  their  physical  shapes.  The  critical 
organs  may  be  represented  by  a few  discrete  points. 

This  approximation  technique  can  simplify  the 
optimisaticn  problem  while  saving  time,  and  can  be  used  to 
develop  a reference  system  to  prepare  standard  guidelines. 
The  method  yields  good  results  which  can  then  be  tailored  to 
the  individual  case. 


After  defining  the  target  volume  and 
treatment  planning  involves  selection  of  t 
that  reduce  the  dose  to  the  normal  tissues 
adequate  target  coverage.  It  is  possible 
beam  parameters  as  the  optimum  variables; 


other  structures, 
s while  providing 


recent  optimitation  methods  use  the  most  sensitive  beam 
parameters  as  the  optimum  variables  to  obtain  favorable 
results  economically.  Many  noncoplanar  arcs  are  used  far 
starsotactic  radiosurgery.  Accordingly,  many  beam 
parameters  (e.g.  the  position  of  Isocenter,  collimator  size 
or  shape,  arc  spacing,  and  the  length  and  weight  for  each 
in  dose  optimization.  The  problems  are 


included 


hov  to  reduce  the  variables  systanatlcally  and  choose  the 
optimum  variables  for  the  particular  situation.  Teats  ware 
done  to  find  the  effects  of  optimum  variables  on  spatial 
dose  distribution  with  noncoplanar  arcs.  The  following  are 
sooe  intereetinq  points  obtained  from  the  tests: 

(1)  The  shapes  of  high  isodose  surfaces  were  almoet 
spherical  with  a single  fixed  isocanter,  even  with  a change 
in  the  position  of  arc  and  weighting. 

[2)  The  isocenter  separation  and  collimator  siae  with 
multiple  isocanters  or  collimator  shaping  were  sensitive 
variables  in  changing  the  high  isodose  shape. 

(3}  The  shape  of  the  beam  profile  or  dose  gradient  was 
dependent  on  the  isocenter  separation  and  collimator  size. 

(4)  Different  are  weighting  for  each  isoeenter  shifted  all 
the  isodose  lines  toward  the  more  weighted  isocenter, 

(5)  The  shape  of  low  isodose  lines  changed  with  a change  in 
arc  direction  and  weighting. 


[6)  The  shape  of  high  isodos 
significantly  with  a change 
position. 

The  tests  stimulated  id 
optimum  parameters.  In  orde 


0 obtain  the  corresponding 
and  collimator  size  must  be 


position 


Ula  active  optUBun  variablaa.  As  an  altamativa,  changing 
the  collimator  shape  could  be  an  affective  parameter  in 
obtaining  the  desired  target  dose  shape.  Arc  variables 
could  also  be  useful  optimum  variables  to  minimize  dose  to 
the  surrounding  normal  tissues  after  a uniform  dose  to  the 
target  volume  is  obtained,  since  those  variables  change  low 
iaodose  lines  without  changing  high  isodose  lines.  In 
Chapter  6 two  different  methods  are  discussed  about 
selecting  beam  parameters  for  obtaining  target  corresponding 
doses  using  a multiple  isocenter  approach  and  a field 
shaping  approach  with  projection  technique. 


3-D  Pose  calculailonal  Model 


After  selection  of  beam  parameters,  dose  distributions 
are  calculated  using  the  exact  geometry  of  the  patient  and 
beam  configurations.  The  3-D  dose  algorithm  is  required  to 
calculate  doses  at  any  point  in  the  patient  for  any 
specified  gantry  angle,  turntable  angle,  collimator  angle 
and  collimator  size.  Since  the  dose  elgcrithm  operates  on  a 
3-D  patient  representation  stored  as  set  of  multiple 

B sectional  contours,  calculating  dose  in  3-D  space 


The  use  of  the  exact  dose  model  in  optimization  is  not 
efficient,  especially  when  ua  utilize  mathematical 
optimization  programming.  A faster  and  more  efficient  dose 


calculation  algorithm  ii 

paramatars  much  faster  \ 
optimisation.  In  Chaptar  5,  a 
Is  davelopad  for  noncoplanar  a 
and  a fitting  tachnlgue. 


luirad  tor  the  optimisation  to 
il.  This  finds  the  optimum 

. a fast  approximate  dose  mode 


Ootinizatton  c 


An  essential  goal  of  radiosurgery  is  to  reach  the 
optimum  dose  distribution  within  ths  target  volume  and 
surrounding  healthy  tissue.  The  most  Important  subject  is 
optimization  criteria  fcr  the  optimal  dose  distribution. 
Criteria  for  optimal  external  beam  radiation  therapy  can  be 
separated  into  two  stages:  clinical  optimization  and 
physical  optimization.  Optimal  treatment  planning  at 
clinical  stages  reguires  a knowledge  of  the  relative 
response  of  tumor  and  normal  tissues  to  irradiation  in  terms 
of  dosage  and  fractionation.  Physical  optimization  consists 
of  the  development  of  a treatment  plan  which  will  produce  a 
preferred  dose  distribution,  our  concern  here  is  physical 
optimization,  since  it  is  difficult  to  obtain  necessary 
biologic  response  Information  on  the  single  fraction 
treatments  uasd  in  radiosurgery  at  present. 

Defining  optimization  criteria  in  a mathematical  form 
and  determining  the  relative  importance  of  different 


criteria  are  the  most  serious  obetacles  to  performing 
autonatefl  optimiaation.  Much  e£  physical  optiffliaation  is 
the  specification  of  naxiisum  and  minimum  dose  to  the  target 
and  maximum  dose  to  normal  structures  at  risk.  Table  1-1  is 
a list  of  some  of  the  quantitative  parameters  that  have  been 
used  to  evaluate  physical  dose  distributions. 


Many  variables  are  used  in  radiosurgery.  It  Is  very 
Important  to  obtain  an  optimal  dose  distribution  within  a 

in  the  first  stage  is  critical  to  the  success  of  real  tine 
optimisation  in  both  visual  and  computer-aided  design 
optimisation. 

The  pre-study  for  the  selection  of  optimum  variables 
must  be  analysed  to  find  the  relationship  between  dose 
distribution  and  optimum  variables.  This  information  must 
be  stored  so  that  it  can  be  used  to  determine  the  optimised 
parameters,  such  as  isocenter  separation  or  collimator  sire, 
to  shape  the  dose  distribution  to  ths  target.  Furthermore, 
these  guidelines  are  useful  in  determining  the  initial 
values  for  the  computer-aided  design  optimisation.  A useful 
relationenip  between  optimum  isocenter  separation  and 
collimator  sire  Is  used  to  fit  an  elongated  target  within 


isodose  surface  in  Chapter  6. 


1-1  ouantttjtlv  Evaluation  at  Pliviol 


- Target  volume  shape 

- so  or  90%  dose  volume 


- Dose  gradient 

- Dose  volume  histogram 

- Treatment  volume 

- Integral  dose 


OPtiaitation  with  eanacralnta  Problen 


The  nain  Ujoitations  In  planning  are  tolerance  doses  of 
neighboring  tissues  and  critical  organs  as  well  as  the 
tolerance  dose  of  the  noraal  connective  tissue  within  the 
tunor  region.  Most  earlier  optinization  methods  were  Based 
on  least  square  fitting  to  sinlnlze  the  variance  of  dose 
distributions  within  the  applied  target  volume,  subject  to 
constraining  doses  at  nearby  critical  structures. 

Eadlosurgeons  may  be  more  interested  in  minimizing  dose 
to  the  critical  organ  with  constraining  points  for  the 
iSdSired  dose  to  the  target,  since  it  is  always  possible  to 
deliver  a lethal  dose  to  the  tumor,  and  a large  amount  of 
dose  is  delivered  to  the  brain  within  a short  time.  An 
additional  constraint  is  the  limit  on  variables,  which  must 
be  specified  with  reason  after  consideration  of  mechanical 
limitations  or  other  parameters. 


Some  optimization  techniques,  such  as  linear  or 
quadratic  programming,  hava  been  applied  to  deal  with 
specific  classes  of  optimization  with  constraints  problems 
in  radiation  therapy.  However,  all  of  these  optimization 
attempts  have  incorporated  a fairly  limited  set  of  linear 
variables  such  as  beam 


weights 


lengthy  recalculation  tine  and  the  difficulty  of  dealing 
with  nonlinear  variables.  Moat  of  the  variables  used  in 

position,  field  size,  and  arc  position,  which  are  very 
ihportant  factors  in  modifying  dose  distributions  to  fit  the 
target  shape  while  avoiding  dose  to  the  critical  organ.  A 
more  general  category  of  algorithms  referred  to  as  nonlinear 
programming  must  ba  considered  for  the  solution  of  the 
optimization  problem  in  LiNAC-based  radiosurgery. 

since  nonlinear  programming  shows  some  disadvantage  in 
time  required  to  find  the  optimum  parameters,  and  it  is 
sometimes  difficult  to  find  the  global  optimum,  using  too 
many  variables  is  not  efficient.  However,  nonlinear 
programming  used  with  a good  estimate  of  initial  variables 
and  a simple  objective  function  could  be  a very  powerful 
technique  in  the  search  for  optimum  variables.  Therefore, 
the  fast  3-0  dose  algorithm  with  a proper  optimization 
criteria,  which  makes  the  objective  function  simple.  Is  an 
important  factor.  The  proper  selection  of  optimum  variables 
which  change  isodose  shape  are  also  important  factors.  Some 
cases  are  discussed  for  obtaining  optimal  dose  distributions 
using  nonlinear  programming  with  constraints  in  chapter  7. 

In  the  following  chapter  we  review  previous  studies  on 
dose  modeling  and  optimization  methods  used  in  radiation 
therapy.  In  Chapter  3 the  system  of  LIHAC-based 
radiosurgery  is  briefly  described.  Dose  modeling  for  single 


and  nultiplB  movin;  beans  is  discuaaad  in  Ctiaptere  4 
Experimental  and  analytical  optimisation  techniques 
discussed  in  Chapters  6 and  7. 


CHhPTEI)  2 

REVIEW  or  LITERATURE 


Photon 


Dose  Model 


Historical  overview 


A central  element  in  any  radiotherapy  treatment 
planning  sysceo  is  an  algorithm  Tor  dose  calculation.  Tsien 
(1955)  is  usually  mentioned  as  Che  TirsC  user  of  computers 


graphically  using  measured  Isodoee  charts,  had  been  the 
standard  method  before  the  possibility  of  utilizing 
computers  vas  studied.  A vide  variety  of  computational 
models  for  the  numerical  evaluation  of  dose  distributions  of 

approaches  taken  by  various  models  can  be  categorized  into 
three  groups:  data  bases  of  measured  values:  empirical  or 
semi-empirical  formulae:  and  analytical  method. 


that  of  determininq,  by  measurement,  the  relative  dose  in  a 
tvo-  or  three-dimensional  array  of  points  In  the  field. 
These  arrays  of  points  and  the  corresponding  relative  dose 
measurements  are  stored  and  the  relative  dose  at  any 
arbitrary  point  is  arrived  at  by  interpolation.  Holmes 


vith  new  e^ipoent,  the  Prograsmed 
: (197XJ  also  measured  beams  for  a 
and  stored  the  m 


(1970)  used  this  me 
Console.  Bentley  t 
limited  number  of  s 

files  on  magnetic  tape  or  disc.  This  method  h 
advantage  of  speed  and  a direct  relationship  b 
measured  and  calculated  doses.  However,  apar 
fact  that  these  methods  require  an  impressive  m 
effort  and  mass  data  storage  facilities,  the  approach  is  not 
applicable  to  the  more  complicated  treatment  beam 
arrangements. 

Because  of  this,  interest  has  grown  in  another  method 
using  empirical  or  semi-empirical  functions.  This  method  is 
based  on  the  application  of  an  empirical  or  semi-empirical 
function  to  describe  the  radiation  field,  sterling  et  al. 
(1964)  developed  this  type  of  model  for  depth  dose  and  cross 


sections. 


were  developed  by  many  other  groups  (Kalnaes  1966:  Orchard 
1964:  Komelson  t Young,  1975:  Thomas  6 Haybittle,  1975), 

The  projective  beam  model,  which  is  more  empirically 
oriented,  has  been  studied  by  the  van  da  Oeijn  group  (Van  de 
Ceijn,  1965,  1986;  Van  de  Galjn  et  al.,  1984).  Another 
popular  model,  where  primary  and  secondary  radiation  are 
calculated  separately,  was  presented  (Worthley  6 Cooper, 
1967:  Gupta  & Cunningham  1966)  and  Introduced  the  concept  of 
scatter-alr-ratio.  The  scatter  air  ratio  concept  was 
extended  in  a study  by  Cunningham  (1972).  A less  amplrieal 


metAod  of  representing  megavoltage  beams  was  developed  by 
Haywood  at  al., (1979) . That  model  paralleled  the 
Cunningham  model  and  took  advantage  of  rectangular  beam 
geometry  to  simplify  calculations.  These  empirical  models 
have  been  used  with  a correction  factor  to  account  tor 
inhOBogeneity.  Before  CT,  the  effective  path  length 
approach  was  widely  used.  The  power  law  tlesue-air-ratio 


method  was  developed  by  Batho  (1964)  and  generaliaed  by 
Sontag  and  Cunningham  (1977),  and  Webb  and  Fox  (1980). 
However, this  method  does  not  explicitly  calculate  the 
ssstter  component.  The  differential  scatter-air-ratio 
method  proposed  by  Cunningham  (1972)  does  explicitly 
separate  primary  and  scatter  but  requires  more  computer 
time.  After  the  development  of  CT,  pixel-based  algorithms 
were  developed  by  several  groups.  These  are  modified  Batho 
power-law  methods  (Cassel  et  al.,  1981 r fuen  at  al.,1984), 
the  equivalent  tissue-air  ratio  method  (Sontag  6 Cunningham, 
1979)  and  the  delta  volume  method  (Hong  et  al,,  1984  and 
Wong  9 Henkelman,  1983).  These  methods  take  advantage  of 
the  three-dimensional  anatomical  information  from  CT  scans 
and  provide  improvements  over  the  earlier  methods. 

Analytical  methods  are  baaed  on  physical  grounds  rather 
than  empirically  chosen  functions.  It  uses  a general  form 
of  the  transport  equation,  with  various  techniques  being 
employed  tor  the  solution.  Boyer  and  Mok  (1984)  attempted 
to  the  computation  of  scatter  dose  by  using 


approach 


Fourier  transfone  techniques  to  convolve  e scatter  dose 
distribution  kernel  with  the  prinarv  fluence  distribution. 
Ahnesjo  (1984)  also  described  an  investigation  of  the 
possibility  of  using  Fourier  transfoms  with  point-spread 
functions  to  calculate  the  dose  distribution  in  external 
photon  bean  therapy.  Mackie  and  Scringer  (1984)  developed  a 
convolution  framework  for  radiation  therapy  dose 
calculations.  The  dose  model  using  Konte  Carlo  photon 
transport  has  been  studied,  with  particular  eaphasia  on  the 
scatter  component  of  such  calculations  by  Hanglos  et  al., 
(1987) . Many  other  analytical  models  are  being  studied  and 
developed,  however  analytical  models  are  found  to  be  less 
useful  for  treatment  planning  because  elaborate  mathematical 

with  the  photoh  beam  model  not  mentioned  here.  This 
literature  was  reviewed  following  the  stream  of  trends. 


contemporary  Bela.tea_s.tudies 


some  experience  with  the  treatment  of  brain  tumors 
using  multiple  moving  beams  shows  the  necessity  and 
efficiency  of  a three-dimensional  radiotherapy  treatment 
planning  system.  Noncoplanar  beans  enter  the  patient  in  an 
arbitrary  configuration  in  stereotactic  radiosurgery.  This 
type  of  treatment  technique  requires  a 3-D  planning  system. 
This  system  must  employ  calculational  techniques  for  four 


individual  model  representations:  patient. 


computational,  treatment  Oeam  setup,  and  visualisation. 

compute  the  3-D  dose  distribution  tor  multiple  beam 
arrangements  by  adding  the  beams  and  combining  the 
relationship  between  the  (our  model  representations. 

Several  dose  algorithms  have  been  described  which  calculate 
three-dimensional  dose  distributions  for  LINAC-based 
radiosurgery  recently  (Hartmann  et  al.,  1984,  1987;  Fibe  et 
al.,  1987).  These  3-D  algorithms  are  based  on  an  isocentric 
model  (Khan  et  al.,  1980),  which  is  one  of  the  empirical 
dose  models  that  describe  the  dose  for  a single  field,  and 
is  used  as  a basis  for  superimposed  fields.  For 
radiosurgery  planning,  an  algorithm  is  needed  that  permits 
the  calculation  of  dose  distribution  (or  any  single  fixed 
circular  field  in  any  orientation,  as  well  as  of  a number  of 
superimposed  fields. 

According  to  Hartman  et  al.  <1984  and  1985],  the  dose 
was  considered  only  near  isocenter.  The  dose  model  of  a 
single  field  was  represented  by  a measured  relative  dose 
file  and  a measured  depth  dose  curve  in  a cylindrical 

approximated  by  that  of  pencil-like  beam.  The  resulting 
dose  distribution  was  calculated  by  superimposing  single 
fields.  Rice  et  al.,  (1987)  used  a simple  isocentric  dose 
model  similar  to  that  used  by  Hartmann  et  al.,  (1985); 


however,  he  used  somewhat  different  measured  faotora.  The 

axis  ratio  were  used  as  basic  beam  data  by  that  group, 
since  those  two  models  actually  ignored  the  inverse  square 
factor,  dose  computation  nay  not  be  accurate  at  points  far 
from  the  isocenter.  Pike  et  al.,  (19B7)  used  a dose  model 
based  on  the  2-D  Milan  and  Bentley  algorithm  (Milan  ( 
Bentley,  1974)  which  uses  measured  central  axis  percentage 
depth  doses  and  beam  profiles  at  several  depths  in  water. 

The  TMR  was  calculated  from  the  percentage  depth  dose  data 
using  Khan's  expression  (Khan,  19S0)  , and  renomalised  to 
the  Isocenter.  This  model  is  a more  accurate  method  of 
obtaining  the  dose  within  the  single  field  than  the  previous 
three  beam  models.  One  disadvantage  of  this  model  is  the 
enormous  effort  for  the  measurements.  Many  other  groups 
have  mentioned  the  use  of  a 3-D  dose  computational  algorithm 
in  their  papers;  however,  they  were  not  described  in  detail. 


The  concept  of  optimiration  as  applied  to  radiotherapy 
planning  is  concerned  with  determining  the  best  possible 
dose  distribution  for  a given  treatment  situation.  Any 
method  which  attempts  such  optimisation  must  incorporate  the 
criteria  which  define  this  ideal  dose  distribution.  A 

out  by 


detailed 


Hope  et  al.,  fl967),  The  usual  starting  point  for 
quantitative  optlaization  waa  the  adoption  of  score 
functions  as  quantitative  measures  of  the  quality  of  the 
treatment  plan.  Hops  et  ai.,  fl967)  introduced  the  ooncept 
of  calculated  Judgement  which  allows  a more  quantitative 
comparison  of  isodose  distributions  to  be  made  through  the 
score  functions.  The  approach  requires  a considerable 
amount  of  data  preparation  and  computer  time  and,  as  it  is 
non-interactive,  it  does  not  allow  rapid  recaloulation  if 
resulting  plans  show  undesirable  distributions. 

Another  approach  using  linear  progranning  was 
investigated  by  several  workers.  Hope  et  al.,  (1972)  and 
Orr  (1972)  attempted  to  optimize  the  physical  or  geometric 
dose  distribution  through  the  use  of  linear  programming, 
which  optimized  scores  that  varied  linearly  from  a desired 
dose  at  a selected  number  of  points.  Bahr  (19S8)  and  Hodes 
(1974)  also  used  linear  programming  to  minimize  the  integral 
dose.  Initial  results  using  this  technique  appear  to  be 
encouraging  but  a problem  does  arise  in  the  choice  of  the 
objective  function  since  this  must  be  linear,  by  definition. 


target  region,  w 


a workers  defined  color 


uniformity  of  dose  to  the  target  region  was  also  defined  by 
setting  upper  and  lower  dose  constraints  at  preselected 
points  within  the  target  area.  Redpath  et  al.,  (1976) 


developed  an  optinization  method  in  which  the  variation  in 
dose  over  the  volume  is  minimized  subject  to  constraining 
doses  at  nearby  critical  structures.  The  calculation 
proceeds  by  way  of  a quadratic  programming  which  minimizes  a 
second  order  objective  function  subject  to  linear 
constraints. 

McDonald  and  Rubin  (1977]  described  an  alternative 
approach  using  a least  square  method,  which  could  be  one  of 
quadratic  programming  similar  to  the  method  used  by  Redpath 
et  al.,  (1976|.  However,  they  used  different  weights  per 
sample  point  to  compensate  for  distortions  due  to  irregular 
spacing  of  sample  points  along  the  target  contour, 
starkschall  (1984J  developed  a constrained  least  square 
optimization  method  incorporating,  as  constraints,  maximum 
allowable  doses  delivered  to  critical  organs,  as  well  as  the 
raquirement  that  all  beam  weights  be  non-negative. 

Those  earlier  methods  all  differ  from  one  another  in 
either  their  choice  of  an  objective  function,  or  their  means 
for  optimizing  the  objective,  or  both.  Those  efforts, 
moreover,  were  limited  to  optimizing  beam  weights  and 
reflected  only  two  dimensional  considerations.  The  field 
sizes  and  bean  orientations,  which  may  require  nonlinear 
programming  (NLP) , are  selected  visually  by  t 
planner,  since  it  is  extremely  difficult  to  quantify  a 
optimization  variables  and  to  find  global  optimum  for 


inothar  approacB  is  visual  optimisation,  whieh  has  been 
the  most  common  method.  The  planner  observes  the  effects  of 
varying  any  of  the  parameters  affecting  the  dose 
distribution  and  so  arrives  Interactively  at  an  acceptable 
distribution.  Houever.  by  a visual  interactive  procedure 
alone,  the  full  effect  of  all  these  parameters  cannot  be 
investigated  because  of  the  time  Involved. 

The  use  of  computer  controlled  systems  with  the  visual 
technique  or,  sometimes,  with  mathematical  programming  were 
studied  to  adjust  field  direction,  width  and  dose  rate. 
Computer-controlled  radiation  therapy,  or  dynamic  treatment, 
is  simply  an  extension  of  the  concept  of  moving-field 
radiation  therapy.  Trumph  et  al.,  (1961)  expanded  this 
technique.  They  added  such  refinements  as  mechanical 
devices  Chat  change  field  size  and  beam  direction  during  the 
treatment,  as  wall  as  synchronous  shaping  and  shielding. 
These  concepts  were  further  developed  by  Takahashi  (1962), 
who  used  a multileaf  collimator  with  independently  movable, 
multiple  pairs  of  jaws  to  define  unsymmetrical , irregular 
fields  (conformation  therapy).  Mantel  et  al.,  (1978) 
investigated  the  possibilities  of  placing  the  therapy 
machine  under  computer  control  during  treatment.  Mantel  at 
al.,  (1978)  discussed  some  ccmputationsl  methods  to 
determine  the  values  of  the  field  size  and  dose  rate  through 
visual  methods  using  linear  programming  as  functions  of  the 
gantry  angle.  As  another  feasibility  study  of  computer 


controlled  redleticn  therapy,  its  application  to  produce  a 
wedge-shaped  dose  distribution  by  noving  collinator  jaws  was 
evaluated.  Here  advanced  cenputer  controlled  radiation 
therapy  technlgues  ware  described  by  Chin  et  al.,  (1961, 
1983).  The  parameters  varied  during  treatment  included  dose 


collinator 


>le  translation, 
They  also  Inplenented  a least  square 


method  Ccr  Che  purpose  ' 


i automatic  optimisation  (Chin  e 


therapy  are  uel 
slow  in  coming, 
complexity  of  t 


f computer  controlled  radiation 

one  reason  for  this  is  the  increased 
e system  and  the  accompanying  increased  ris)c 


of  error.  However,  computer  controlled  radiation  therapy 
using  multileaf  collimators  has  been  developed  for 
optimization  of  more  general  moving  beam  treatments  (Talcai 
et  al.,  1987;  Lind  & Brahme,  1987;  Bower  et  al.,  1985). 


In  recent  years,  advances  in  computer  and  graphics 
display  technologies  have  opened  a vide  range  of 
possibilities  for  3-D  treatment  planning.  Conformal  therapy 
using  multilaaf  and  beam's  eye  view  display  techniques 
(frass  et  al.,  1967;  Goitein  et  al.,  1983;  Mohan  at  al.. 


1988)  is  most  sffactivB  to  shape  target  volume  in  3-0  space, 
as  well  as  in  3-0  compensator  design,  while  analytical 
optimization  approaches  have  been  applied  in  teletherapy  or 


brachytherapy  optiBization  using 


ling  (TaJcai 


at  al.,  1987)  or  quadratic  programming  (Toll!  fi  Ragnhult, 
1987;  chung-bin  b Malyala,  1987;  Tabush  at  aX.,  1988). 

several  preliminary  studies  have  been  propoaed  in 
recent  proceeding  papers,  and  research  is  ongoing  to  solve 
more  complicated  optimization  problems.  Dynamic  treatment 
(varying  treatment  parameters)  (Schlegal  at  al.,  1987),  or 
treatment  under  temporal  variation  oC  one  or  more  treatment 
parameters  during  the  irradiation  is  one  of  the  approaches. 

A multiaxial  treatment  technique  was  suggested  by  Aohi  et 
al.,  (1987)  in  order  to  overcome  the  disadvantage  of 
conformal  treatment  planning  with  a multileaf  collimator. 

The  technique  is  based  on  a summation  of  partially  rotated 
conformal  techniques  with  different  iso-axes  for  multileaf 
multiple  planning  objects  to  which  patient  volume  is 
divided.  Hosoi  et  al.,  (1987)  described  a dose 

volumes  were  put  in  one  original  target  volume,  and  the  dose 
distributions  from  each  planning  object  were  calculated  with 
conformal  rotation  techniques.  A dose  distribution  nearly 
corresponding  to  the  shape  of  the  oriqinal  target  volume  was 


At  present  no  detailed  method  is  mentioned  in  the 
literature  for  the  optimization  of  dose  distribution  for 
stereotactic  radiosurgery. 


DESCRIPTION 


CHAPTER  3 

OF  STEREOTACTIC  PHOCEDDRES 

Stereotactic  radiosurgery,  in  various  forms,  is  in  use 
in  numerous  countries  (e.g.  France,  Italy,  Sweden,  United 
States,  USSR,  West  Germany)  for  the  treatment  of  a variety 
of  lesions.  It's  utilisation  and  applications  are  steadily 
increasing.  In  this  chapter,  the  procedure  of  a 
stereotactic  radiosurgery  system  la  described  briefly, 
including  unique  mechanical  and  software  features  of  the 
University  of  Florida  radiosurgieal  system  (Friedman  c Bova, 

In  order  to  explain  the  stereotactic  procedure,  the 
three  steps  of  the  procedure  (localisation,  dose  planning, 
and  treatment)  must  be  examined  separately.  The  ultimate 
accuracy  of  the  full  procedure  is  dependant  on  each  of  these 
steps  and  on  the  consistency  of  the  approach. 


A Iccallzation  system  is  an  essential  prerequisite  for 
stereotactic  radiosurgery.  The  localization  technique 
Includes  the  stereotactic  localization  of  the  target  in  the 
brain  by  using  a head  frame,  and  determination  of  the 


ooordinatas  of  the  target  point  with  reference  to  this  frame 
and  adjustment  of  the  target  point  to  the  isocenter.  This 
is  accomplished  by  one  of  two  means.  Currently,  the  method 
of  choice  is  through  stereotactic  angiography.  The 
procedure  begins  with  the  stereotactic  ring  being  fitted  to 
the  patient.  A specially  designed  intracranial  localizer 
box  (or  angiographic  localising  device)  then  rigidly 
attaches  to  the  Brown-Robert-Bell  (BHH)  Head  Ring.  This 
device  has  bean  described  in  other  literature  (Bidden  S 
Barth,  19B7)  and  consists  of  four  sets  of  reference 
alignment  markets.  The  coordinates  of  the  sixteen  markere 
are  precisely  known  with  respect  to  the  BRW  coordinate 
system  (or  stereotactic  ring) . Two  localization  films  are 
required  for  reconatruotlon  of  the  target  position  and 
determination  of  the  target.  An  analysis  of  the  projected 
target  and  radio-opaque  oar)ters  permits  reconstruction  of 
the  target  position  and  determination  of  magnification, 
while  this  part  of  the  procedure  allows  the  coordinates  of 
the  target  relative  to  the  laoeliaation  ring  to  be 
determined,  more  anatomical  information  is  needed  for 
dosimetry  analysis. 

The  next  step  replaces  the  angiographic  localizing 
device  with  another  localizer  specially  designed  for 
localization  in  CT.  This  is  the  standard  BRW  CT  Localizer. 
The  BRW  head  ring,  with  CT  localizer  attached,  is  clamped  to 


parallel 


head  ring.  Five  nillimeter  coneecutive  scans  are  obtained 
troB  above  to  below  the  head,  If  we  can  sea  the  target  in 
CT.  Three  nillineter  scans  are  obtained  from  above  to  below 
the  target  while  S nb  scans  being  dona  for  renaining  part  of 


the  head,  if  we  can  not  see  the  target  in  CT.  If  the  target 
volume  can  be  identified  in  the  CI  image,  then  the  xys 
coordinates  of  the  target  volume  are  again  calculated.  A 
simple  algorithm  that  transfcriiis  spatial  points  in  CT 
coordinates  to  BRW  stereotactic  frame  cccrdinatea  has  been 
presented  in  a recent  paper  (Saw  et  al.,  1987).  with  the 
digitally  encoded  data  from  the  CT  scan  and  the  two 
angiographic  films,  the  data  is  then  transferred  to  the 


dosimetry  computer  system.  The  CT  scan  provides  three- 
dimensional  anatOBical  information  of  the  patient,  allowing 
a solid  patient  model  to  be  constructed.  The  coordinates  of 
the  target  volume  from  the  angiogram  and  the  CT  scan  data 


The  next  concern  is  a dosimetry  system  to  evaluate  the 
dose  delivered  to  the  target  and  normal  tissues  in  the 
patient,  vrtiile  it  generates  an  optimal  dose  distrlbutlcn 
that  will  satisfy  a set  of  clinical  criteria  for  the 
patient.  Radiosurgery  treatment  in  the  brain  requires 
detailed  information  on  three-dimensional  dose  distribution. 


t planning  program  is  a 
prerequisite  for  treatment  plan  optimisation.  It  must  cover 
3-D  methods  for  representing  the  patient,  the  dose 
distributions,  and  the  beam  settings.  The  method  is  baaed 
upon  a patient  anatomical  representation  which  is  obtained 
from  CT  images.  The  ideal  system  for  administering 
radiation  therapy  is  to  deliver  a tumorioidal  dose  to  the 
target  volume,  not  to  irradiate  tissues  outside  that  volume, 
and  to  do  this  with  relative  ease  and  efficiency.  Using  the 
3-D  computation  algorithm  and  an  appropriate  optimisation 
method  (visual  or  analytical),  optimum  beam  setting 
parameters  can  be  obtained.  The  target  position,  spacing 
between  arcs,  the  size  of  the  collimator,  and  the  variation 
in  arc  length  and  weight  will  be  parameters  used  in  the 
optimization.  The  procedure  needed  for  dose  planning  in 
stereotactic  radiosurgery  is  shown  in  figure  3-1. 


Radiation  Treatment 


After  dose  planning,  the  radiation  treatment  system  i 
assembled  and  tested  after  the  photographic  method  of  Lutz 
et  al.,  (198B).  A S KV  isocentric  LIHAC  with  a 100  on 
source-to-axis  distance  is  used  at  the  University  of 
Florida.  The  required  additions  or  modificationa  to  the 
LIKAC  to  make  it  useful  for  radiosurgery  are  extra 
colliioator  systems  (to  define  the  small 


circular  fields) 


Figure  3-1.  The  procedure  for  oprinal  dose  planning  in 
stereotacCic  radioaurgery. 


mschanical  accuracy 


a three-baaring  systam  (to  improve  the 
of  the  LIKAC)  (Friedman  & Bova,  1989) . 


One  bearing  controls  the  isocentric  accuracy  of  the 
collimator.  The  second  bearing  controls  the  rotation  of  the 
head  of  the  BRW  floorstand.  These  two  bearing  systems  alone 
cannot,  however,  produce  the  desired  accuracy  if  rigidly 
attached  to  the  LINAC  head.  A gimble-type  bearing  with  a 
sliding  collimator  mount  was  developed  to  avoid  any  torque 
transfer  from  the  LINAC  head.  This  allows  Che  radiosurgery 
mechanical  system,  and  not  the  LINAC  head,  to  determine  the 
isocentric  treatment  accuracy  (Figure  3-2) , 


A series  of  cerrobend  collimators,  IS  cm  in  length  and 
ranging  in  diameter  from  1 to  3 cm  in  steps  of  2 mm 
(including  special  sizes  of  5 mm  and  3.S  cm)  were 
constructed.  The  collimators  have  multi-step  cylindrical 
apertures  which  define,  at  the  machine  isocenter,  circular 
fields  with  the  diameters  constructed.  The  multi-step 
openings  are  tapered  to  match  beam  divergence  and  minimize 
penumbra.  The  accelerator  collimator  is  sat  to  a field 


The  Mchanieal  system  Illustrated  with  lateral 
(upper)  and  overhead  (lower)  (From  the  published 
paper  by  Freedman  and  Bova,  1989).  The  LINAC 
head  moves  around  the  head.  The  patient  is  moved 
to  a new  couch  position  after  each  individual  arc 
IS  completed.  The  rotation  of  the  floor  stand 
: jantry  (B)  are  controlled 

gimbal  bearing  ( 


During  treatneni  the  patient's  head  is  positioned  and 
innobilised  independently  of  the  radiotherapy  couch  using 
the  8RH  floor  stand  (Keilbrun  et  al.,  1983).  The  stand's 
three  orthogonal  drives  and  verniers  are  used  to  position 


CHAPTER  4 

DOSE  HODSL  FOR  A FIXED  SINGLE  BEAM 

Before  we  discuss  the  two  main  approaches  for 
optimlzaticn  of  LINAC-baaed  radiosurgery,  we  must  present 
some  information  on  dose  oodelllng.  In  this  chapter  we 
discuss  the  dose  model  and  beam  measurement  for  a fixed 
single  beam.  This  is  basic  background  for  developing  an 
efficient  3-D  dose  algorithm  for  multiple  moving  beams. 

The  3-D  dose  algorithm  for  multiple  moving  beams  will  be 
considered  in  the  next  chapter. 

A simple  empirical  isocentric  model  (van  Fanthalean, 
19B6I  Rice  et  al . , 1987;  Khan  et  al.,  1980)  will  be  modified 
to  describe  the  dose  for  both  circular  and  rectangular 


circular  Field 


An  isocentric  model  for  single  circular  treatment  is 


BOF<C)  X THB(w,d)  X (SAD/STD)' 


[4 


n point  of  intscest  in  a msdium 

c B collimator  size  defined  at  SAD 
STD  = source  to  target  distance 

u - field  size  at  point  of  interest  m 
expressed  by  w - C (STD/SAD) 

d - depth  of  point  of  interest  m 


r ■ off-axis  distance 

SAD  » source  to  axis  distance  ■ 100  cm 


D^f  ■ the  dose  for  the  reference  set-up 
ROF  =■  relative  output  factor  defined  by 

D(C,  STD-100,  d,,r»0)/D(C,^,SSD“100,d..O) 
TMR  - tissue  maximum  ratio  defined  by 
0(w,d)/D(w,d,l 

OAR  - off-axis  ratio  defined  by 

0(C,STD,d,r)/D(C,STD,d,r-0) 


The  measurement  data  for  tissue-maximum-ratio  (THR) 
off-axis-ratio  (OAR)  are  the  main  function  data  for  the 
model  developed.  These  measurements  are  identified  as 
beam  data.  Early  work  on  basic  bean  data  included  two 
parts.  The  first  part  was  to  determine  the  detector 


system  foe  our  smsll  bean  measurement.  The  second  part  vae 
to  eomulate  basic  beam  data  Iron  our  limited  measurements. 

While  we  tested  many  detectors  (e.^.,  iso-octane  liquid 
ion  chamber,  Markus  ion  chamber,  TLO,  film,  and  diode),  the 
diode  detector  was  chosen  as  our  primary  beam  measurement 
device  due  to  the  qood  characteristics  and  the  small  size  of 
the  sensitive  volume  (2  mm  diameter).  Flqure  4-1  shows  good 
detector  response  while  representing  the  linear  relationship 
between  diode  reading  and  dose  for  various  x-ray  beam 
energies.  Figure  4-2  shows  the  good  agreement  between  the 
diode  and  other  detectors  (fils  and  TLO)  for  the  OAB 

The  effort  to  formulate  the  percentage  depth  dose, 
output  factor,  and  beam  profile  was  made  for  the  normal 
rectangular  collimator  size  in  the  previous  reports  (Dixon 
et  al.,  1977;  van  Dyk,  1977;  (ioudek  et  al.,  1963;  Johns  i 
Cunningham,  1983).  The  formulation  of  the  THR  for  the  small 
circular  collimator  was  described  in  a recent  paper  (Rice  et 
al.,  1967).  In  order  to  formulate  both  the  TMR  and  OAR,  a 
limited  number  of  dose  measurements  were  obtained  from  the 
reference  set-up.  Then  the  possible  empirical  formula  were 
used  to  fit  the  data.  Those  formulae  were  later  extended 
and  generalized  for  any  real  treatment  setup  by  finding  the 
relationship  between  fitting  parameters  and  beam  parameters. 
Empirical  formulae  for  6 MVX  at  Radiation  Therapy  at 
Gainesville  Incorporated  (RTGI)  and  8 MVX  at  Shands  Teaching 


Detector  response 
energies. 


number 


off-axis  ratio  for  field  sizes  of  1.  2.  and  3 
cm  (SAD*100  cm)  scanned  at  depths  of  l.S  (aj, 

10  (b),  and  5 (c)  cm  in  solid  water  using  film, 
diode,  and  TLD. 


Hospital  were  derived  to  represent  THR  and  OAR  for  any  real 
treatment  setup  fron  a limited  number  of  dose  measurements. 


The  TKR  was  measured  fron  the  fixed  source  to  chamber 
(or  source  to  target)  distance,  lOQ  cn  with  the  diode  and 
build-up  of  solid  water.  Figure  4-3  shows  trends  of  TMR 
data  measured  (data  points]  for  two  different  x-ray  beams. 
TMR  data  were  fit  to  a function  of  the  font  beyond  the 
build-up  region  (d^-1.5  cn  for  6 HVX  or  2 cm  for  8 MVX) : 

TKR(STD,C,d)  - A(C)  X exp(-«(C)  x d)  (4-2) 

where  u(C]  is  the  effective  linear  attenuation  coefficient 
in  water  for  collinator  field  size  C,  and  A(C)  allows  one  to 
account  for  the  shift  in  d^  with  field  size.  The  curves  in 
figure  4-3  represent  the  fit  of  TMR  data  tor  1,  2,  3 and  4 
cm  collimator  sizes.  Figure  4-4  shows  a plot  of  A and  s 
values  as  a function  of  collimator  field  size.  The  values 
for  A and  p are  scattered  in  a linear  fashion.  Figure  4-5 
shows  plots  of  TMR  vs.  depth  tor  different  source-to-target 
distances.  These  show  a good  agreement  (to  within  1%) 
between  the  TMR  obtained  from  the  different  source-to-target 
(STD)  distances,  since  TMR  is  independent  of  source  to 
target  distance  and  dependent  on  field  size  w^  at  depth  d, 
the  formula  for  TMR  can  be 


expressed  by 


Ti.ssue*naxinuiD-ratio  for  6 HeV  x-ray  (a)  and  s 
HeV  x-ray  (b) . The  points  denote  the  neaeured 
bean  data  and  the  curves  denote  the  analytic 
function. 


6 MV  X-ray 


6 MV  X-ray 


8 MV  X-ray 


Tissue  naximuio  ratio  for  different  source-to- 
distance  (STD) . The  points  represent 


THR(w.d) 


A(W)  - A,  t A,  w 

«(w)  - + .1,  w 


w ■ C X STD  / SAD,  SAI>-100  cm 


[4-3) 

(4-4) 


(4-S) 

(4-6) 


and  A,,  Aj,  (j,  and  (ij  are  coeffielente  at  linear  rectrassion. 

Figure  4-6  shows  the  TMR  trends  at  or  near  d^  for  6 
HVX  for  both  STD-100  and  STI-90  using 
values  of  TKR  for  the  larger  field  si: 
than  those  for  the  smaller  field  size 
d^,  and  vice  versa  beyond  d^.  it  is 


filB  measurament.  The 
e are  slightly  smaller 
in  the  region  below 
possible  to  derive  an 


empirical  fonnula 
the  difference  of 
is  small,  it  is  desirable 
for  all  the  collimators  in  the  build-up  region.  Figure  4--; 
shows  the  THR  data  and  cubic  polynomial  fit  within  the 
build-up  region. 


TMR  within  the  build-up  region.  Since 
among  the  different  collimator  sizes 
obtain  an  approximate  formula 


Shends  Teaching  Hospital 


gives  trends  similar  to  those  of  6 
at  Gainesville  Incorporated  (RTGl) . 


Table  4-1  Table  Ot  Fittina 
Ilaaue  Maximum 


Oft-  Axis  Ratio 


001  ■ axs 


Ud-up  reqion  for  t 
axtterent  energies,  s and  a MaV  x-ray  using 
diode.  The  points  denote  measured  data,  and 
curve  denote  cubic  polinooial  fit. 


oeaaured  with  the  diode  at  d^  and  at  lo  co 
r and  6 MVX  beans  while  using  varying 

;e-to-target  distances. 


depth  for  6 MVX  i 
combinations  of  collimators  and  i 
Figures  4-e  show  the  result  of  01 
different  collimator  field  sires 


since  the  difference  between  OAAs  at  depth  at  any  cff- 
axis  distance  is  small,  any  profile  dependence  on  depth  can 
be  ignored  for  normal  clinical  applications.  The  curves  in 
figures  4-8  represent  multi-depth  fits  for  each  collimator 
using  the  modified  Cunningham  model: 


o«t(STD,c,r)  - 1 - 0.5  X exp[-o,  x (wyj  - r)/p 

-B  X (wyi  - r)VP'] 

=■  t ♦ (0.5  - t)  X axp(-«,x  (r  - w/3)/p] 

for  r > wy2  (4-7b) 

where  w„  is  given  by  Eg.  4-6,  and  r represents  off-axis 
distance.  Penumbra,  p,  is  expressed  by 


(4-8) 


collimator  distance  (70 


Off-axis-ratio  for  6 Mev  X-ray  (a)  and  8 HaV  X- 
ray  (b) . The  points  represent  the  measured 
beam  data  for  collimator  sizes  of  1,  2,  and  3 
cm  at  depths  of  d^  and  10  cm  with  fixed  STD 
of  100  cm.  The  curves  denote  fitting  formula. 


6 MV  X-ray 


8 MV  X-ray 


Figure  4-9  shows  OAR  beam  data  at  and  10  c 
with  the  respective  formula  generated  curves  far  t^ 
different  STDs  with  a collimator  size  of  2 cm.  Th« 
accuracy  of  fit  was  obtained  from  fitting  formula  i. 
collimators,  since  transmission  does  not  change  much,  the 
constant  value  of  0.02  was  used  for  t.  Most  of  the  fitting 
curves  agree  well  with  measured  data  except  for  smaller 
collimators.  The  addition  of  a guadratio  term  to  the 
exponent  for  the  inner  penumbra  description  gives  an 
improvement  to  field  sizes  smaller  than  1.4  cm.  This  makes 
it  possible  to  create  profiles  which  have  a flat  maximum  in 
the  center,  while  maintaining  t 
the  penianbra  region  for  the  sma 
Figure  4-10  shows  trends  of 
and  Og  with  the  variation  of  STD 
different  beam  energies.  There  i 
determine  fitting  parameters  for 
STD.  One  direct  method  i 
parameters.  Another  is  t 


e appropriate  gradient  in 
1 oollimator. 
the  fitting  parameters  o. 
nd  collimator  size  for  t 
a several  options  to 
ny  collimator  size  and 
o interpolate  between  known 
e relationship  between 


fitting  parameters  and  beam  parameters  such  as  collimator 
size  or  STD.  Flgurs  4-11  shows  composite  linear  regression 
of  a,  and  Cj  with  STD  for  multi-collimators  ignoring  the 
collimator  dependence  of  a,  and  Cj.  Figure  4-12  shows 
anothsr  composite  linear  regression  for  multi-collimatore 
for  s MVx  with  film  measurement.  The  formula  for  OAR  now 
can  ba  generalized  with  new  expressions  for  a.  and  a : 


ray  (b] . The  points  represent  the  eeasured 
data  at  d_  and  10  cn  depth  for  three 
different  STDs  with  a collimator  size  of  2 cm. 
The  curves  denote  fitting  formula. 


6 MV  X-ray 


8 MV  X-Ray 


Fitting  parameters  o£  OAH  function,  and  e^ 
va  source-to-target  distance  for  different 
collimator  settings.  (a)  for  6 MeV  X-ray, 
(b)  a,  for  6 MeV  X-ray,  (c)  a,  for  6 MeV  X- 
ray,  (d)  a,  for  6 MeV  X-ray. 


6 MV  X-ray 


8 MV  X-ray 


Fitting  parameters  and  a,  vitb  composite 
linear  regression.  The  thin  solid  lines 
represent  the  trends  of  Individual  fitting 
parameters,  a,  and  a,  for  different 
collijaator  sires  and  STDs.  The  thiclc  solid 
lines  denote  the  composite  linear  regression 
for  multi -collimators . 


8 MV  X-ray 


4-12.  Pitting  paramatBrs  a,  and  Oj  vs. 
target  diatance  for  colliaator 

X-ray,  (b)  Oj  for  6 Mev  x-ray. 


Sourca  Tsfjei  Dtslanee  (Ctrl 

(a) 


Q,  ♦ Qj 


(4-9) 


n additional  expression  for  small  collimators  (C  < 


Q,  * Q,  STD 


The  coefficients  of  regression  are  shown  In  table  4-1.  A 
formula  may  be  derived  from  the  biquadratic  form  which 
includes  dependence  on  both  collimator  size  and  STO. 
However,  this  is  not  recommended  as  more  parameters  only 
provide  a minimal  Increase  in  the  accuracy  of  the 
expression. 


Relative  Output-Factor 

To  express  the  relationship  between  dose  and  monitor 
units  for  different  collimator  sizes,  relative  output 
factors  (ROF)  are  used.  The  relative  output  factor  was 
measured  at  d^  with  source-to-chamber  distance  of  100  cm 
with  the  reference  setup  (source  to  surface  distance  of  100 
cm  at  d^  for  field  size  of  10  X 10  cm).  Figure  4-13  shows 
the  measured  ROF  for  different  collimator  sizes.  The  ROF 
data  nay  be  fit  to  a function  of  the  form: 


(4-11) 


Collimator  Field  Size  (mm) 


Relative  output  (actors  vs.  eolZiioator 
settings.  Collinator  settings  define  the 
field  size  at  isocentec.  The  points  denote 
neasured  data,  the  curves  denote  least  square 


single  circular 


D(0,STD,d,r)  - D,.,  X 80F(C)  X TMH(w,d)  X (SAD/STD)' 

X OAR(C,STD,r)  (4-12) 

where  all  function  values  are  explicitly  represented  by  the 
formulae  developed  here. 


Rectangular  Field 


An  isocentric  model  for  a single  rectangular  treatment 
is  given  by 

D.(H,,W2,STD,d,y;,t,)  - 0,,,  X ROF(H,,Hj)  X TM8(W,,Wj,d) 

X (SAD/STD)"  X OAS(W,,Wj,SrD,d,yj,Zj)  (4-lS) 

where  W.,,  and  are  field  sites  measured  at  isocenter,  and 
Wf  and  Wj  are  field  sizes  at  the  point  of  interest.  The 
off-axis  distance  is  represented  by  two  off-axis  variables, 
y^  and  z^,  in  the  cartesian  collimator  coordinate  system, 
other  parameters  are  the  same  as  those  for  the  circular 
bean,  with  circular  field  size  replaced  by  rectangular  field 
size.  The  rectangular  field  is  used  for  Che  feasibility 
study  of  optimal  dose  distribution  associated  with  conformal 
optimization  using  beam's  eye  view  in  chapter  6.  The 


rectangular  (laid.  Instead  of  saking  enomous  efforts  to 
obtain  the  basic  bean  data  (or  rectangular  fields,  we  used 
basic  beam  data  for  circular  fields  in  representing  basic 
beam  data  for  rectangular  fields  using  equivalent  field 


Tlssue-Maximum-Batle 

From  the  relationship  between  square  fields  and 
circular  fields  {Sjamgard  t slddon,  1962) , the  equivalent 
circular  field  diameter,  w,  can  be  represented  by 


where  w,  is  a size  of  the  perimeter  of  the  square  field, 
kny  rectangular  field  size  can  be  converted  into  an 
equivalent  square  field  by  using  a conversion  expression 
(Johns  & Cunningham,  1983) : 


From  Eq.  4-14  and  4-15,  the  equivalent  circular  diameter,  v, 
is  expressed  in  terms  of  rectangular  field  size: 


'niR  for  rectangular  fields  can  be  obtained  from  TKR  for 
circular  fields: 

T«R(w,,Wj,d)  - TKB(w,d)  (4-17) 

where  w is  given  by  Eq.  4-lS,  and  the  analytic  tom  of 
TKR(Vjd)  is  given  by  Eq.  4-3  to  4-6. 

Off-Axls-Batio 


The  OAB  can  be  represented  by  two  separate  off-axis 
ratios  as  in  a primary  beam  profile  (Johns  & Cunningham, 

OAR(W,,Wj,STQ,d.yj,*j) 

- OAK(H,,STD,d,y^f  X OAR(Wj,STD,d,  (4-18) 

For  the  OAR  for  each  axis  of  rectangular  bean,  the  OAR 
expression  developed  for  circular  fields  may  be  used 
approximately,  ignoring  depth  dependence  as  well  as  for 
circular  fields: 

OA»(W,,STD,yj)  - OAR(C,STD,r)  (4-19a) 

OAR(Wj,STO,*j)  “ OAR(C,STD,r)  (4-19bl 


given  by 


single  rectangular  fields  is 


0.(W,,Wj,STO,d,y,,s,)  ■ D,,,  x HOF{W,,W,)  x TMH(w,,Wj,d) 

X (SJU3/STD)’  X OAa{W,.STD,yj)  x OAH[Hj,STD,  z^)  (4- 

where  all  function  values  are  given  by  Egs.  4-16,  17,  19 
and  19b,  and  are  explicitly  represented  by  the  fomula 
developed  for  circular  fields. 


MODEL 


CHAPTER  5 

TKRBE-DIHERSIOHAL  DOSE 

MULTIPLE  HOVIHS  BEAMS 

In  chapter  4 we  discussed  dose  models  for  fixed  single 
beams.  Multiple  noncoplanar  arcs  are  used  for  LIHAC-based 
stereotactic  radiosurgery.  Exact  3-D  dose  models  for 
non-coplanar  arcs  were  developed  using  fixed  single  beam 
models  for  both  circular  fields  and  rectangular  fields. 

Those  exact  3-D  dose  models  will  be  the  bases  to  evaluate 
optimal  dose  distributions  through  experimental  optimization 
techniques  such  as  the  multiple  isocenter  technique 
(circular  fields)  or  field  shaping  (rectangular  fields) . 

Past  approximate  dose  models  were  also  developed  using  the 
exact  dose  model  for  both  coplanar  and  multiple  noncoplanar 
arcs.  These  fast  3-D  dose  models  will  be  used  to  find  the 
optimum  beam  parameters  through  computer  based  analytical 
optimization. 


Dose  Model 


Circular  ?leia 

slnola  treacment  b«an  aodel  for  angle  orientation 

From  the  single  isocentric  dose  model  given  by  Eg. 
4-12,  the  calculation  of  the  dose  at  each  point  of  interest 
D requires  the  depth  d,  the  off-axis  distance  r,  the  source 
to  target  distance  5TD,  and  the  collimator  siae  C with 
gantry  and  turntable  orientation  pj.  From  these 
parameters  the  field  size  v at  each  point  is  easily 
calculated  by  the  relationship:  w - c (STD/SAD).  When  we 
use  different  isocenter  positions  for  different  arcs,  the 
depth  d and  off-axis  distance  r of  each  point  on  the  dose 
grid  for  each  increment  of  gantry  rotation  8,  in  the 
algorithm  are  the  only  parameters  which  need  to  be  derived 
from  the  patient  contour  and  isocenter  position  Ij  for  each 
arc  j . The  formula  to  express  the  dose  at  defined  point  m 
with  two  additional  variables  6|,  dj  can  be  written  by 

D(C|,STD,j,,d|j,,r|jJ  = D,,,  X BOF(C,)  X TMR(w,,., d,,,) 

X (SAO/STO|,,) ' X OAB(C,,STD|j„r,.,)  (S-1) 

e,  » gantry  angle  orientation 

dj  = turntable  angle  orientation 


as  in  Eg.  4-1. 


other  parameters 


defined 


The  beam  setbing  paranebera  musb  be  debermined  first  to 
calculate  the  dose  at  the  paint  of  interest.  They  consist 
of  collimator  else  C,  burn  table  orientation  gantry 
angle  orientation  9|,  the  frame  coordinates  of  the  isccenter 
I.  for  each  arc  j,  and  treatment  time  t.  Using  the  3-D 
computation  algorithm  and  optimisation  method,  optimum  beam 
setting  parameters  can  be  obtained. 

Dose  grid 

For  dose  calculations  a matrix  of  points  is  established 
within  the  image  being  viewed  through  the  display  plane.  A 
cartesian  3-D  dose  grid  is  the  most  popular,  however  a dose 
matrix  on  an  oblique  line  or  oblique  plane  may  be  desired  In 
certain  cases.  Three  cases  of  dose  grids  are  discussed: 
cartesian  reference  recti-linear  BRW  coordinates  relative  to 
the  BRH  head  ring  (axial,  sagittal,  coronal  plane  etc.): 
oblique  line  dose  grids  defined  by  a line  passing  through 
two  arbitrary  points  A and  B,  or  defined  by  angle  parameters 
a and  pj  and  an  oblique  plane  defined  by  three  arbitrary 
points  A,  B and  C,  or  defined  by  a rotation  of  a cartesian 
reference  coordinate  system  about  a cutting  point  M by  angle 
parameters,  a and  0. 

Xt  the  dose  matrix  and  beam  setting  parameters  are 
determined,  the  calculation  of  the  dose  on  the  dose  grid 


requires  d,j,,  r,j,,  STO,j,,  and  from  our  dose  nodel  in  Eg. 
5-1.  The  method  is  based  primarily  on  the  definitions  of 
the  coordinate  systems  and  transformations,  or  vector 
analysis  in  three-dinensional  apace. 

Pimad  machine  coordinate  avstaa.  The  fixed  coordinate 
system  is  an  orthogonal  coordinate  system  vhich  is  fixed  in 
the  treatment  room  (Figure  5-1) . The  origin  is  at  the 
isocenter.  The  x,  axis  is  directed  vertically  upwards,  and 
the  Zf  axis  is  horizontal  and  directed  towards  the  gantry. 
The  gantry  rotates  about  the  z,  axis  in  x,-y,  plane.  The 
turntable  rotates  about  the  x^  axis  in  a plane  parallel  to 
the  y,-z,  plane.  The  angles  and  d,  define  the  positive 
directions  of  rotation  of  the  gantry  and  turntable.  The 

This  coordinate  system  will  be  the  basio  coordinate  system 
for  finding  the  un)cnown  parameters. 

Beam  coordinate  system.  The  gantry  bean  coordinate 
system  is  a cylindrical  coordinate  system  (Z,r,d)  which 
rotates  about  the  Zf  axis.  The  c*0  axis  corresponds  to  the 
beam  axis  and  the  origin  is  the  intersection  of  the  beam 
axis  and  the  isocenter.  The  * coordinate  is  actually  not 
required  since  we  can  assume  symmetry  about  * direction. 

Isocentric  BRH  coordinate  system.  The  isocentric  BRH 
coordinate  system  (Xp,yp,Zp)  is  defined  by  a rotation  of  the 
Fixed  machine  coordinate  system  (x^.y^.z^)  by  di  about  the  Zf 
axis  in  the  fixed  coordinate  system.  The  origin  is 


75 

isocenter.  During  rotation  of  the  gantry  and  turntable,  the 
patient  is  at  rest  with  respect  to  the  isccentric  BRW 
coordinate  system. 

BBW  coordinate  system.  The  BRW  coordinate  system 
(x,,y,,2,)  is  the  stereotactic  frame  coordinate  system.  The 
isocenter  positions  and  dose  grid  for  display  plane)  are 
based  on  this  BRW  coordinate  system.  The  BRW  coordinate 
system  is  the  reference  coordinate  system  to  represent  the 
individual  patient  geometry.  The  CT  slices  are  transferred 
to  the  BRW  coordinate  system. 

The  coordinate  transformation  from  [VJ^  to  [V),  is 
given  by 

£V]p  - CV],  - [«|],  (S-2a) 

where  X|  is  the  translation  vector  which  translates  [V]|  to 
[V]  , and  represents  the  vector  position  of  the  isocenter,  I 
in  the  fV]g  system.  Eg.  5*2a  can  be  written  as 


* yi  - Vi 


(5-2c) 


The  coordinate  transfomation  from  the  isocentric  BRW 
coordinate  ayaten  to  the  fixed  coordinate  systen  is  given  by 

[V],  - R,(«,)  [V]„  {S-3a> 

where  is  the  rotation  operator  which  rotates  the  [V} 

system  by  the  rotation  of  turntable  about  Xp  or  Xp.  Using 
the  explicit  expression  of  rotation  operator  (Siddon,  19S1) , 


The  unknown  parameters  must  be  determined  on  a defined 
dose  grid  for  any  isocenter  position  and  angle  orientation. 
We  will  consider  three  cases  of  dose  grids,  cartesian  1<*D 
BRH  coordinates;  oblique  lines;  and  oblique  planes. 

on  cartesian  BRM  coordinate  dose  grid.  From 
geometrical  relations  in  Figure  5-2  between  a 
line  IjMj  (I,K,  or  IjMj) , 


Figure  5-2. 


coordinate  systems  vith^’cartesian  dose  grids 
(dot  line)  in  recti-linear  BRW  coordinate 
system. 


X, 


X 


>Ss: 


(5-4) 


where  l^Hj  and  I’jk  can  be  expressed  in  herns  of  vector  units 
in  the  fixed  coordinate  systen: 


I,.M  - X,  i * y,  j + *,  k (5-6) 

Rearranging  Eg.  5*>4  with  Eqs.  5-5  and  5-6: 


Eg.  5-7  can  be  expressed  in  terns  of  the  isocentric  BRH 
coordinate  systen  with  Eg.  5-3: 

[5-8) 

Eg.  5-8  can  be  expressed  in  terns  of  the  BBH  coordinate 
systen  with  Eg.  5-2: 


where  (X.,yj,zp  are  the  coordinates  of  the  isocenter  Ij  for 
arc  j in  the  BRW  coordinate  system. 


Eq.  S-9  is  a generalized  equation  to  detemine  Z,j,  in 
the  cartesian  BRW  coordinate  s/steia.  The  dose  grid  m can  be 
defined  in  any  cartesian  BRH  plane  (e.g.  Xg-y|  for  the  axial 
plane,  x,-z,  for  the  sagittal  plane,  yi-z^  coronal 

plane,  etc.].  If  is  found,  un3cnown  parameters  such  as 
d,  r,  STD,  and  v can  be  determined  from  their  relationship. 

The  depth,  d,.,  is  determined  by 


= p.jH  = d|,|  - Z,j.  (5-10) 

where  d|j|  is  the  depth  of  isocenter,  and  Z(j^  is  given  by  Eq. 

The  problem  of  finding  the  depth,  d|j,  is  now  reduced  to 
finding  the  depth  of  isocenter,  d,j,  at  each  increment  of  the 
gantry  rotation  e,  and  turntable  rotation  about  isocenter 
position  Ij.  The  depth  of  isocenter  is  dependent  on  the 
individual  patient  contour,  isocenter  position,  and 
transformation  of  CT  data  to  BRW  coordinate  system,  etc. 

To  calculate  the  depth  of  the  isocenter,  d^j,,  a simple 
algorithm  has  been  developed  which  is  based  on  a database 
containing  a sampling  of  points  on  the  surface  of  the 
patient’s  head  using  a ray  tracing  technique  (Siddon,  1985] 
by  the  UP  research  team  and  other  institutes  (Pike  et  al., 
19B7) . 

However,  for  the  research  purpose  of  optimization,  a 
spherical  head  model  will  be  used  as  our 


reference  patient 


nodel,  and  nhe  coordinate  of  origin  of  the  sphere  will  be 
assumed  to  be  (0,0,0)  in  the  BRW  coordinate  system.  Then, 
the  depth  of  isocenter  will  be  found  for  a given  gantry 
angle,  turntable  angle  and  isocenter  position  by  their 
mathematical  relationship.  The  point  P on  the  surface  of 
the  spherical  head  shown  in  figure  5-3  can  be  expressed  in 
terms  of  the  BRW  coordinate  system  given  by 


X,  - R cos  e,  + X,  (5-11) 

y,  - R sin  e,  cos  o,  + Yj  (S-12) 


where  R is  the  depth  of  isocenter  for  a given  gantry  angle 
Si,  turntable  angle  and  isocenter  position  (Xj,Y,,Zj)  in 
BRH  coordinates.  The  square  sun  of  Eqs.  5-11  through  5-13 
corresponds  to  the  square  of  radius  of  the  sphere  R,  and  is 
given  by 

X,'  * y,'  + z,‘  . R„- 


» 2|  sin  e,  sin  *.)  ♦ (X,-  + Yj'  + 2,-)  (5-14) 

If  we  arrange  Eq.(S-14),  we  obtain  a quadratic  equation: 


(5-lS) 


(5-16) 


c . (X,‘  * * Zj')  - V (5-17) 

The  solution  of  this  quadrstic  equation  nay  be  written  in 


where  b and  C are  given  by  Eq.  5-16  and  S-17.  From  the 
geoaetrical  relationship,  the  off-axis  distance.  r,j,  is 
determined  by 


expressed  with  some  algebra  by 


*^i|»  ■ C V * V * 'b'  ■ *!]■'  (5-20) 

Eq.  5-20  can  be  expressed  in  terms  of  the  BRW  cartesian 
coordinates,  with  Eq.  5-2  as: 

■ H’S  - ^i)'  ty.  * y,)'  + (2.  - - Zi|.‘  ]’'* 

(5-21) 

where  Z,.^  is  given  by  Eg.  5-9.  The  source  to  target 
STD,j^  is  determined  by 


distance, 


detenfilned  by 


I (STO,j,  / SAD) 


where  c.  is  the  collimator  field  sise  defined  a 

Bqs.  5-10,  16,  17,  18,  21,  22,  and  23  now  constitute 
an  explicit  set  of  algebraic  equations  for  the  determination 
of  unknown  paraneters  on  the  cartesian  BRW  coordinate  dose 
grid  with  our  spherical  model  (radius  Rq)  in  terms  of  the 
known  parameters  d,,  0j,  c,  and 

on  oblique  line  dose  arid.  He  may  be  interested  in  the 
dose  distribution  on  any  oblique  line  through  any  two 
arbitrary  points  A and  B,  especially  when  we  consider  using 
two  isocenters  to  treat  the  inclined  elongated  target. 

We  define  a dose  grid  on  the  oblique  line  AB 
(rig.S-4a).  The  first  point  A will  be  the  origin,  and  the 
second  point  B will  determine  the  direction  of  the  oblique 
line.  When  A and  B are  determined  in  BRW  coordinates,  the 
position  vector  k on  the  line  AB  can  be  expressed  in  terms 
of  two  position  vectors  )£,  and 


(S-24) 


Figtire  5-4.  Geonetrical  relationship  between  BPW  coordinate 
systejo  and  oblique  line  dose  grid;  (a)  defined 
by  a line  passing  through  two  arbitrary  points 
A and  B,  (h)  defined  by  a origin  point  A and 
angle  parameters  a and  B. 


in  BRW  coordinates  (5-2S) 

* C^/yb'^bl  coordinates  (5-26] 


s - (x,.y,.*,)  <s-27) 

where  t is  a scalar  value  representing  the  distance  fron  the 
origin  A to  point  X along  the  line  AB.  We  define  the  scaled 
dose  grid  t'  as: 


ilfe  - X.I  - J(x„  - X,)-  ♦ (y^  - y,)-  ♦ (i„  - t.)‘ 

The  position  x defined  in  cartesian  BRW  coordinates  is  now 
redefined  in  terms  of  the  scaled  dose  grid  t'  on  the  oblique 
line  AB  using  Bqs.  S-24  to  5*29. 

X,  - X.  + t'  (Xj  - X.)/|Xb-25.l  (5-30) 

y,  - y,  + t'  (yj  - y.)/lx,-J£.l  (5-31) 

2,  =.  X.  + t'  (2„  - 2, 


(S-32) 


The  oblique  line  can  be  also  defined  by  an  oblique  angle  o 
and  from  the  axis  passing  through  the  origin  A 
(Fig.S-4b).  Then,  the  position  s in  the  BSW  coordinate 
system  is  expressed  by 


(5-33) 


(S-3S) 


where  t"  is  another  dose  grid  on  the  oblique  line  defined  as 
above.  £qs.  5-33  to  5-35  are  useful  to  represent  the  radial 
dose  distribution  in  any  radial  direction  defined  by  a and  0 


used  as  an  explicit  set  of  algebraic  equations  for  the 
determination  of  the  unknown  parameters  on  the  oblique  line 
dose  grid,  which  is  defined  by  Eqs.  5-30  to  S-32  or  by  Eqs. 
5-33  to  5-35  with  our  spherical  model  (radius  Rg)  in  terms 
of  the  known  parameters  9,,  C,  and  Ij. 

On  oblique -Plane  dose  arid.  Ue  may  be  interested  in 
dose  distributions  on  any  oblique  plane  through  three 
arbitrary  points  A,  B,  and  C,  especially  when  we  consider 
more  than  two  isocenters  to  treat  an  inclined  irregularly 
shaped  target.  This  is  convenient  when  we  consider  a plane 


containing  both  the  target  area  and  critical  organ  of 

when  we  consider  an  oblique  plane  defined  by  three 
points  A,  B and  C in  figure  S-5,  the  position  X on  that 
plane  can  be  represented  in  two  ways;  cartesian  reference 
BRH  coordinates,  or  oblique  plane  coordinates  defined  by 
scalars  s and  t.  Suppose  three  points  A,  B and  C are 
deter&ined  in  the  reference  BRW  coordinate  system  and 
represented  by  vectors  a,  b and  e; 

« ’ (x,,y,,z,)  (S-36) 

6 - (Xj,y^,Zj)  (S-37) 

= - (x.,y.,t.)  (S-38) 

A dose  computation  algorithm  was  developed  to  calculate 
dose  in  the  reference  cartesian  coordinate  system.  With 
this,  oblique  plane  coordinates  should  be  expressed  in  terms 
of  cartesian  BRW  coordinates  in  order  to  find  unknown  bean 
parameters  on  the  oblique  plane  defined  by  the  three  points, 
using  vector  notation  the  position  X oan  be  expressed  in  the 
cartesian  BRW  coordinate  system  by 

X = OX  = X,  i t y,  j -HZ,  k (5-39) 

or  can  be  expressed  in  the  oblique  plane  coordinates  defined 
by  two  orthogonal  axes  s and  t by 


Geometrical  relationship  between  BRH  coordinate 
system  and  oblique  plana  dose  grid  defined  by 
three  arbitrary  points  A,  B,  and  C. 


uhsre  i,  } and  k are  unit  vectors  of  three  orthogonal  axes 
in  the  cartesian  BRH  coordinate  system,  and  s and  t are  unit 
vectors  of  two  orthogonal  axes  In  the  oblique  plane 
coordinate  systen. 

The  vectors  s and  t are  expressed  by 


(5-41) 


OQ  - OP 
lOQ  - 0P| 


Iq  - Pl 


a expressed  by 


(5-42) 


(5-43) 


The  vector  p can  be  expressed  by 


(5-45) 


(5-46) 


Sinplifying  the 
coordinates: 


of  cartesian 


where  a and  b are  given  in  Egs.  s-36  and  37,  and  |b  - a|  is 
represented  by 


l*>  ' a|  = + (y,-y.)‘  + (a„-z.)‘  (S-48) 


represented  by 


J + (y,-y,)'  * (t.-z,)- 


(S-50) 


q = OQ  - OA  + AQ 


(5-43a) 


given  by  Eg. 


given  by 


where  e is  given  by  Eg. 


is  given  by 


09  - OP  <i 

log  - 0P|  I? 

where  q and  p are  given 
is  represented  by 


(5-43a) 


and  |q  - p| 


k - Pl  ■ * tVP»>'  * f"!.-?.)'  (S-5i) 

The  conponents  at  all  vectors  are  now  expressed  in 
tens  of  cartesian  BRW  coordinates.  Fron  Eqs.  5*39  and  S-40 
the  position  X expressed  by  the  s-t  coordinate  systesi  on  the 
oblique  plane  defined  by  three  points  A,  B and  C is  now 
redefined  in  terns  of  the  cartesian  SAW  coordinate  systen 
= 

X,  - X.  -I-  8 s.  + t t.  (5-52) 

y,  = y,  ♦ s 8,  t t t,  (5-53) 

3,  = 2,  + s s,  + t t,  (S-54) 

now  will  be  used  as  an  explicit  set  of  algebraic  equations 


for  the  determination  of  unknown  beam  parameters  (depth  of 
isocenter  Zj.,.  depth  of  point  of  interest  d,,  , off-axis 
distance  r|j,,  source  to  target  distance  STDu,  and  field  size 
at  point  of  interest  w,.,)  on  the  oblique  plane  dose  grid 
(s-t)  defined  by  three  points  A,  B and  c fron  known 
parameters  (isocenter  position,  collimator  size,  turntable 
rotation  and  gantry  rotation) . 

However,  a more  convenient  way  to  represent  the 
direction  of  the  oblique  plane  might  be  to  use  two  angle 
parameters,  a and  B,  as  described  in  the  oblique  line  dose 
grid.  The  oblique  plane  dose  grid  is  defined  in  the  rotated 
BRW  coordinate  system  [V),,  which  is  defined  by  a rotation 
of  the  BHW  coordinate  system  by  a about  any  x-  line  passing 
the  origin  point  A to  an  intermediate  system,  followed  by  B 
about  its  own  y axis  (Fig. 5-6). 

First,  we  define  the  translated  BRW  coordinate  system, 
(V),  fron  the  BRW  coordinate  system,  [VJ,  with  translation 
vector  s,  (or  OA) : 

(V)t  “ [V],  - [3S.),  (5-55) 

The  T and  B systems  are  parallel.  The  origin  of  the  T 
system  is  A,  ie.  (X^,V^,Z^)  in  the  BRW  coordinate  system. 

Second,  we  define  a rotating  BRW  coordinate  system,  [V], 
with  the  rotation  angle 


[V],  coordinate  system 


Geonetrical  relationship  between  four 
coordinate  eystens  with  oblige  plane  dose  grid 
defined  by  a cutting  point  A and  angle 
parameters  a and  5. 


paranaters  a and  8-  Using  rotation  operators, 


transformation  from  [V],  to  [V],  is  therefore 

(V],  - H,[a)  [VJ,  (S-S6) 

Eg.  5*56  can  be  also  expressed  using  the  unitary  operator 
property  by 


CVIt 


B,(-a)  [V], 


(5-57) 


where  R«(e]  and  tty(8)  are  given  by 


(S-5B) 


R,(fl)  - I"?  1 V"  *• 

'sin  8 0 cos  8\  (5-59) 

The  dose  grid  is  now  defined  on  the  x,-y,  plane  in  the 
(V),  system,  which  is  rotated  from  the  x,-y,  plane  in  the 
(V),  system  by  a about  the  X,  axis  and  8 about  the  y,  axis. 
The  position  x defined  in  cartesian  BRW  coordinates  is  now 
redefined  in  terms  of  the  dose  grid  x,-y„  on  the  rotated 
oblique  plane  with  Eqs.  5-55  to  5-59! 

[V3,  - R,(-Q)  IV(-fl)  [V),  + [a.],  (5-60) 


Simplifying 


matrix  algebra. 


(S-61> 


y,  . y,  cos  a * X,  sin  o sin  fl  + Y.  {5-s?) 

z,  > y,  sin  a + x,  cos  a sin  i9  ■»  Z.  (S-G3) 

Eqs.  5-9,  10,  16,  17,  18,  21,  22,  and  23  now  will  bs 
used  as  an  explicit  set  of  algebraic  equations  far  the 
determination  of  unknown  parameters  on  the  oblique  plane 
dose  grid,  which  is  defined  by  ggs.  5-si  to  5-63  with  our 
spherical  model  (radius  Rq)  in  terms  of  the  known  parameters 


The  final  algorithm  to  express  dose  at  the  point  of 
interest  n in  the  defined  dose  grid  for  a noncoplanar  arc  j 
which  may  have  different  isccenter  positions  is  expressed  in 
the  form  given  by 


where  is  given  by  Eg.  (5-2) , 

The  unknown  beam  parameters  for  one  gantry  (9,)  and 
turntable  orientation  (pj)  of  each  arc  j are  determined  in 
the  cartesian  reference  BhW  coordinate  system,  or  the 
oblique  plane  coordinate  system.  TKR  and  OAR  values  are 
then  obtained  by  empirical  fitting  formulae  developed  in 
Egs.  4-2  through  4-11  with  the  beam  parameters  determined 


orientation 


Proa  the  single  isocentric  dose  model  given  by  Eg. 
4-13,  the  calculation  of  the  dose  at  each  point  of  interest 
a requires  the  depth  d,  the  off-axis  distances  and  z^, 
the  source  to  target  distance  STD,  the  collimator  size  U, 
and  Wj,  and  field  size  v and  vith  gantry,  turntable,  and 
collimator  orientation  9,,  0j  and  The  field  size  and 

at  each  point  is  easily  calculated  by  the  relationship, 
u,  = (STD/SAD).  When  we  use  different  isocenter 
positions  for  different  arcs,  the  depth  d and  off-axis 
distance  y^  and  z^  of  each  point  on  the  dose  grid  for  each 
increment  of  gantry  rotation  6,  in  the  algorithm  are  the 
only  parameters  which  need  to  be  derived  from  the  patient 
contour  and  isocenter  position  X-  for  each  arc  The 
formula  to  express  the  dose  at  defined  point  m with  three 
additional  variables  e^,  o*,  and  can  be  written  by 

D,„  X BOF{W„j, X TMH{w„,i.,,W„,i.,,d,j.) 

X (SAD/STD, „)'X  OAR(«,„,,Wj,i„STD„„y.,„>„,Z,„,.^,) 

(S-65) 


0f  » gantry  angle  orientation 

dj  « turntable  angle  orientation 

■ collimator  angle  orientation 


and  the  other  parameters  are  defined  as  in  the  rectangular 
field  in  Chapter  4. 

Beam  setting  parameters 

The  beam  setting  parameters  consist  of  collimator  size 
table  orientation  gantry  angle 
orientation  collimator  angle  orientation  r^,  the  frame 
coordinates  of  the  isccenter  I.  for  each  arc  j , and  arc 
weighting. 

Dose  grid 

Three  casee  of  dose  grids  can  be  used  as  defined  in  the 
circular  field:  cartesian  reference  BRH  coordinates;  oblique 
line  dose  grids;  and  oblique  plane  dose  grids. 

The  calculation  of  the  dose  at  the  dose  grid  requires 

in  Eg.  5-65.  four  coordinate  systems  were  considered  to 
determine  unknown  beam  parameters  for  the  circular  field. 

The  cylindrical  beam  coordinate  system  with  the  rectangular 
field  can  be  replaced  by  the  cartesian  gantry  coordinate 
system  and  collimator  system  (Figure  5-7] . 

The  gantry  system  [V]j  is  defined  by  a rotation  of  the 
fixed  machine  coordinate  system  by  e,  about  the  axis. 

The  origin  is  at  the  isocenter.  The  transformation  from 
[V],  to  [V]j  is  given  by 


[V],  - R,(a,)  [V], 


(5-66) 


The  collimator  system  [V]^  is  defined  by  a rotation  of 
the  gantry  coordinate  system  by  about  the  x,  axis.  The 
origin  is  at  the  isocenter.  The  treatment  field  size  H,  and 
W,  are  measured  at  isocanter.  The  transformation  from  (V), 
to  [V)j  is  given  by 

CV]c  - (Vlj  (5-67) 

The  isocentric  BRH  coordinate  system  [V],  is  defined  by 
a rotation  of  (V),  by  -o*  about  the  x,  axis.  The  origin  is 
at  the  isocenter.  The  transformation  from  (V],  to  (V)^  is 
given  by 

[V),  - R.(-d,)  [V],  (S-68) 

The  reference  BRW  coordinate  system  C^]e  defined  by 
a translation  of  the  [V]|,  system  by  vector  1^: 

tV],  =•  tv]p  + tlj3.  [S-69) 

From  Egs.  5-66  to  5-69,  (V],  can  then  be  transformed  to  the 
[V]p  system  by  the  combined  transformation: 


[V)c  = R,(r, 


.(ej  ( [V],  - tij], 


(5-70) 


(S-71) 


-h;  %i\ 


t fy,  - yj) 


(5-75) 


«c  - (*.  - Xi)  =08  a,  + 


with  the  spherical  head  nodel,  depth  of  isocenter,  d,jj  is 
given  by  Egs.  5-16  to  5-lS: 


,/(b'  - o)  for  H > 0 


(5-18) 


(5-16) 

c - (X,'  + y,'  Zj')  - R,<  (5-17) 

and  y^  and  represent  off-axis  distance  from  the  canter  of 
the  bean: 


y,  = f-(x,  - Xj)  sin  a,  * { (y,  - V,)  cos  <t,  ^ 

(z,  - Z.)  sin  »j)  cos  a,]  cos  r,  + (-{y,  - Y,) 
sin  •,  ♦ (z,-  Z,)  cos  ) sin  t,  (S-76) 

*e  ■ - t-(*«  - >«))  8i'<  8|  + ( (Yg  - Y))  =08  + 

(2,  - Zj)  sin  d,)  COS  a,]  sin  r,  t {-(y,  - Y.) 
sin  d,  + (z,-  Z,)  COS  0,  ) cos  i,  (5-77) 

The  source  to  target  distance  is  detencined  by 

STD|j.  - SAD  - X,  (5-78) 


where  SAD  is  equal  to  100  cm. 

The  field  size,  w and  w^  are  determined  by 


(STO,„ 


(S-79) 


/ SAD) 


(S-BO) 


where  H,  and  are  collinator  field  eizes  defined  at 
isocenter. 

Eq.  5-16  to  S-ia  and  5-74  to  5-80  now  constitute  an 
explicit  set  of  algebraic  equations  for  the  detemination  of 
un]aiown  parameters  on  the  cartesian  BRH  coordinate  dose  grid 
with  the  spherical  model  (radius  Rg)  in  terms  of  the  Icnovn 
parameters  S,,  r,,  w,,  and  1^. 


Using  beam  data  and  the  exact  dose  model,  a study  of 
the  spatial  dose  distribution  for  ccplanar  and  noncoplanar 
arcs  was  carried  out.  The  dose  distribution  generated  by 
the  exact  dose  model  could  be  represented  by  a simple 
approximate  analytic  form  which  is  convenient  and  very 
efficient  for  calculating  dose  distribution  iteratively  in 
the  optimization  procedure.  Two  Xinds  of  simple  dose  models 
were  developed  for  small  collimators  to  represent  the  dose 
in  3-D  space  in  the  patient  coordinate  system,  one  is  a 
spherical  dose  model  for  four  noncoplanar  arcs  about  a 
single  isocenter  with  equal  arc  spacing  (e.g.  the  standard 
arc  used  in  stereotactic  radiosurgery) , and  the  other  is  a 
cylindrical  dcse  model  for  full  or  partial  rotational  arcs. 


Fast  Aoproximate 


HQdel 


spherical  Dose  Model 


When  we  use  multiple  symiDetrical  arcs  with  a single 
isocenter,  dose  distributions  are  approxlnately  spherically 
symnetrical  about  the  isocenter  with  sone  differences  due  to 
the  different  angle  position  (Fig.  5-a) . Our  approach  to 
finding  a closed  form  result  for  approximately  representing 
the  radial  dose  distribution  is  to  utilize  the  nonlinear 
least  square  fitting  (HLLS)  technique.  Among  the 
competitive  analytic  forms,  one  simple  modified  gaussian 
form  was  chosen  for  small  collimators,  which  are  useful  for 
multiple  isocenter  treatments.  The  analytic  form  for 
multiple  arcs  with  equal  spacing,  which  are  commonly  used  in 
stereotactic  radiosurgery  is  given  by 


where  r is  radial  distance  from 
collimator  diameter  (less  than  1.6  cm).  The  p parameters 
are  fitting  parameters.  The  parameters  obtained  by  WLLS  for 
smaller  collimators  are  given  in  Table  S-1.  The  fits  are 
shown  by  the  thlc)c  lines  in  Figure  5-9.  The  thin  lines  are 
radial  dose  distributions  for  different  angle  orientations, 
and  are  determined  by  the  exact  dose  model.  The  gaussian 
form  does  not  fit  well  with  large  collimator  sizes.  The 
more  general  form  was  found  to  represent  radial  dose 
distribution  for  any  collimator  size.  Two  separate  analytic 


directions^for  four  full  rotational  arcs  (a) 
and  four  100'  partial  rotational  arcs  (b) , with 
equal  arc  spacinq  4S".  The  numbers  shown  in 
figure  represent  the  angle  of  radial  direction 
from  zero  gantry  position. 


Fuil  Rotational  Aro 


Partial  Rotational  Arc 


Table  5-1  Paramatarg  for  Approxlaatg  Doae  Model 


in  figure  represent  the  angle  of  ra<3ial 
direction  from  taro  gantry  position. 


Full  Rotational  Arc 


Partial  Rotational  Arc 


forms  are  used  to  represent  the  redial  dose  distribution  for 
both  the  inner  and  outer  penumbra  regions,  and  are  given  by 


Dose,  - 1 - s,  exp[-  s,  x (c/s  - r)  - s,  x (c/j  - r)’  ] 
for  r S C/2  (5-02) 

- s,  ♦ (1  - s,  - s.)  exp(-  s,  X (r  - c/2)] 

The  s parameters  are  fitting  parameters.  The  parameters 
obtained  by  NLLS  for  collimator  sizes  are  given  in  Table  s- 
1.  The  fits  are  shown  by  the  thick  lines  in  Figure  5-10. 

The  thin  lines  represent  the  steepest  and  shallowest  radial 
dose  distributions,  and  are  determined  by  the  exact  dose 
model.  From  examining  the  influence  of  the  simple  spherical 
dose  model,  it  is  concluded  that  it  is  important  to  use 
approximate  forms  with  reasonable  limits  (probably  down  to 
20  t Isodose).  Since  high  dose  is  required  on  the  target, 
the  spherical  dose  model  is  useful  to  shape  the  field  to  the 
target  using  the  multiple  isocenter  approach. 

For  dose  calculations  a matrix  point  is  established  in 
cartesian  BRH  coordinates.  As  an  example,  using  a 1 cm 
collimator,  a dose  model  for  a single  isocenter  with  four 
standard  aros  is  given  by 


r/(0.. 


C))  (5-8fl) 


Figure  5-10.  Radial  dose  profile  for  various  collimator 
size  (1,  2 and  3 om) . The  dark  solid  line 
denote  analytic  fit.  The  thin  solid  line  and 
dashed  line  represent  the  exact  dose  data  of 
the  steepest  and  shallowest. 


where  radial  distance  r is  expressed  in  terms  of  the 
reference  BHW  coordinate  system  (rf.  exact  dose  model)  with 
isocenter  position  I. 

r = [(X,  - X,)'  * [y,  - y,)‘  + (z,  - z, )']''*  (5-as) 

A dose  distribution  for  multiple  isocenters  is  obtained  by 
siopiy  adding  each  contribution  from  each  isocenter.  A dose 
model  for  multiple  isocenters  is  expressed  in  the  form  given 
by 

9 ■ 2 e,  ,5-86) 


where  Dj  is  given  by  Eqs.  5-S4  and  5-86.  Figure  S-ll  gives 
a montage  of  the  spatial  contours  of  dose  distribution  (on 
the  coronal  plane  in  cartesian  BRW  coordinates)  from  single, 
two  and  three  isocenters  with  a l cm  collimator.  The  beam 
profile  AB  and  CD  in  Figure  5-llb  and  c represents  a good 
result  with  a beam  profile  obtained  from  the  exact  dose 
model  (Fig.  5-12). 

Cylindrical  Dose  Model 

The  spherical  dose  model  may  not  be  vary  useful  when 
we  consider  the  optimization  problem  with  constraints  to 
normal  structures.  The  spherical  dose  model  assumes 
symmetrically  spaced  standard  arcs  without  allowing  the 
individual  movement  of  each  arc.  A more 


general  approach 


Figure  5-11.  laodose  distributions  calculated  by  a 

spherical  dose  nodal  on  the  plane  through  the 
isocenters;  single  (a),  two  (b^c),  and  three 
(d) . The  isodoae  lines  displayed  are  fron  20 
to  90  in  10  increments  and  normalized  by 
maxinun. 


Dose  profile  through  A-B  (a)  and  C-D  (b)  in 
Fig.  5-11.  The  doae  profiles  were  calcula 
by  a exact  dose  model. 


ited 


(a) 


Verncai  Axis  (cmi 

(b) 


must  allow  the  movement  of  individual  arcs  and  different  arc 
weighting  to  avoid  significant  dose  to  the  critical  organs, 
which  may  be  located  away  from  the  target  area,  Brahme  et 
al.,  (1962)  developed  an  analytical  approach  for  a solution 
of  an  integral  equation  encountered  in  rotation  therapy.  In 
the  present  study  we  develop  an  empirical  cylindrical  dose 
model  which  may  be  used  to  represent  the  dose  distribution 
for  one  full  rotation  arc  and  one  partial  rotational  arc. 
Since  the  cylindrical  dose  model  does  not  ignore  an  angle 
dependence,  that  model  is  more  accurate  than  the  spherical 
dose  model  over  the  entire  region  including  both  the  target 
area  and  normal  structure  area  of  interest. 

Full  rotational  .are 

One  simple  analytical  form  was  found,  which  well 
represents  both  radial  and  axial  dose  distribution  for  one 
full  rotational  arc  with  small  collimators.  After  a tedious 
search  for  fits  to  a collection  of  200  points  of  exact  dose 
data  we  found  good  results  with  the  equation: 

D,  - q,  X exp((-q2  X * q,  «|Zl-  q.  X r)/  C] 

♦ q,  X axp((-q,  x f - q,  x t‘)/  C ] (S-87) 

where  r and  2 are  radial  and  axial  distances  from  isocsnter 
in  cylindrical  geometry  and  C is  collimator  diameter  (lesa 
than  1.6  cm)  . Dj  is  the  dose  contributed  from  arc  j.  The  g 
parameters  are  fitting  parameters.  The  parameters  obtained 


by  NLLS  for  small  colliaabor  size  are  given  in  Table  5*1. 
Figure  5-13  represents  a aeries  of  radial  dose  distributions 
for  different  axial  positions  with  an  exact  dose  model 
(solid  curves)  and  with  the  cylindrical  dose  model  (dashed 
curves)  using  analytic  fits.  However,  this  continuous 
analytic  form  does  not  fit  veil  with  large  collimator  sizes. 

The  radial  distance  r and  axial  distance  Z in  the 
cylindrical  arc  coordinate  system  are  expressed  in  terms  of 
a fixed  machine  coordinate  system  (rf.  Figure  5-14). 

f = Cx,'  + y,'  )''*  (5-88) 

Z - z,  (5-89) 

The  coordinate,  which  represents  angle  orientation  in  the 
cylindrical  coordinate  system  is  not  required  since  radial 
distribution  shows  symmetry  about  the  angle  direction. 

TO  calculate  dose  in  the  patient  coordinate  (i.e.  BRW 
coordinate)  system  the  dose  matrix  is  transformed  from  the 
Fixed  coordinate  to  the  Patient  coordinate  system.  As 
described  in  the  exact  dose  model,  the  coordinate 
transformation  from  the  isocentric  BRW  coordinate  system  to 
the  fixed  coordinate  system  is  given  by 

[V],  - R.(dj)  [VI, 


(S-: 


X2J. 


Radial  dosa  profiles  for  various  axial  olanes 
paraUel  to  centrall  arc  plane  with  1 on 
collinator  size.  The  dashed  curves  represent 
dose  data  with  an  exact  nodal  and  solid  curves 
denote  cylindrical  dose  nodal. 


coordinate  systen  {x,,  y,,  a,)  centered  ar  the 
isocenter  to  the  cylindrical  coordinate  system 
(r,  ♦,  Z). 


where  ieocentric 


coordinates 


expressed 


terms  of  the  reference  BRW  coordinate  system. 

[VJp  - :V],  - [X|],  (5_2a, 

The  transformation  from  (V),  to  (V),  is  therefore 

[V],  - R.(*j)  ( [V],  - ) (S-90) 

A dose  distribution  for  multiple  arcs  is  obtained  simply  by 
adding  each  contribution  from  each  arc.  A dose  model  for 
multiple  arcs  is  expressed  in  the  form  given  by 


where  Dj  is  given  by  Eqs.  5-87.  Figure  S-15  gives  a montage 
of  the  spatial  contours  of  the  dose  distribution  (on  the 
coronal  plane  in  cartesian  BRW  coordinates)  from  two  arcs 
with  isocenter  position  (0,0,0)  and  collimator  size  1 cm. 

The  isodose  shape  in  Figure  S-15a  shows  good  agreement  with 
that  obtained  from  the  exact  dose  model  (Figure  5-15b) . 

Figure  5-16a  gives  a montage  of  spatial  contours  of  dose 
distribution  on  the  major  arc  plane  from  one  full  rotation 
arc  with  isocenter  position  (0,0,0)  and  collimator  size  1 
cm.  The  isodose  shape  in  figure  S-iea  shows  good  agreement 
with  that  obtained  from  the  exact  dose  model  (Figure  S-16b). 


ircs'^(full^ 

1 by  cylindrica 


In  the  prevlaus  section  a simple  approximate  form  was 

arc  is  used,  the  isodose  curves  tend  to  be  circles  on  the 
arc  plane  (Flg.5-17a).  If  partial  rotation  is  used,  this 
symmetry  is  lost  and  the  region  of  highest  dose  is  shifted 
from  the  axis  of  rotation  away  from  the  missing  sector 
{Pig.S-l7b  to  17f) . If  the  missing  sector  is  a small 
fraction  of  the  360'  rotation,  the  effect  is  negligible.  As 
the  sector  becomes  larger  the  distortion  also  becomes 
larger.  This  1s  true  no  matter  what  the  shape  of  the 

He  adapted  the  cylindrical  dose  model  for  full 
rotational  arcs  to  the  development  of  a more  useful  dose 
model  for  a partial  rotational  arcs.  He  accomplished  this 
by  transforming  the  cylindrical  dose  model  for  full  rotation 
arcs  to  the  modified  cylindrical  dose  model  for  partial 
rotation  arcs  with  an  arc  correction  factor,  which  is 
dependent  on  varying  the  length  and  position  of  the  arc. 

When  we  consider  partial  rotation  arcs,  we  obtain  an 
asymmetry  factor  by  comparing  the  dose  distribution  obtained 
from  the  cylindrical  dose  model  for  full  rotation  with  that 
obtained  from  the  exact  dose  model  tor  partial  rotation 
arcs.  The  arc  correction  factor,  g,  is  now  defined  by 


Dp 


(5-9S) 


Figure  5-17.  Isodose  distributions  on  the  arc  plane, 

calculated  by  exact  dose  nodel  with  various 
partial  rotational  arcs.  (a)  full  rotational, 
(b)  25-,  (o)  SO-,  (d)  100-,  (e)  ISO",  and  (f) 
260- . The  isodose  lines  displayed  are  froa  90 
to  10  in  10  decreaents  and  nomallzed  by 
naximua  dose. 


f 


full 


Eqs.  5-87.  Here,  we  consiaer  a 1 cn  collinator  for  the 
deDonstratlon.  Dp  is  the  dose  from  the  exact  dose  node!  for 
partial  rotation  at  the  same  position.  Since  it  is  not 
practical  to  find  a correction  factor  at  every  position,  It 
Is  desirable  to  find  the  approxinate  correction  factor  over 
the  major  arc  plane  for  the  possible  partial  rotational  arcs 
{100’,  ISO',  etc.).  The  arc  correction  factors  for  50', 
lOQ'  and  180’  partial  rotation  arcs  were  found  for  the 
different  angle  positions  about  five  isodose  lines  (20,  00, 
50,  70  and  80)  of  full  rotational  arcs  (dashed  circles  in 
Fig.  5-18a  to  c) . 

for  100'  and  ISO'  partial  rotation  arcs  are  well  represented 
by  a simple  quadratic  form 

g “ 90  9i  b t gj  t*‘  (S-93) 


ji  = cos  0 (5-94) 

where  e is  defined  as  the  angle  from  the  position  of  the 


133 


S-19.  Arc  correction  factors  tor  five  isodose^lines, 
arcs  SO'  (a),  100'  (b) , and  180'  (c) . 


The  conversion  from  the  fixed  machine 
coordinate  system  (x,,  y,,  z,)  to  the  bisector 
coordinate  system  (x,,,  y,,,  z,,). 


Fitting  parame 
g,,  g,,  and  gj 

rotational  arc 


B of  arc 


function 

full 


139 


100  Arc  Length 


180  Arc  Length 


netsiive  Oose  i%) 


Table  5-2  Coaffleii 


ot  Lii 


Raareaaien  for  g factor 


9l  - 9ii  + 0, 

- 9!i  + 9jj  D, 


%i  ^02  3ii  9i2  9ti 

-0.327  1.663  0.114  -0.116  2.725 


(S-98) 


(5-99) 


To  calculate  dose  in  the  reference  patient  coordinate 
systes  (i.e.  BRH  coordinates),  the  dose  aatrix  is 
transformed  from  the  Bisector  coordinate  system  CV]^,  to  the 
Fixed  Machine  coordinate  system  [V]^,  and  then,  from  (V],  to 
the  Isocentric  BHW  coordinate  system  [V]  , and  finally,  from 
[V]  to  the  reference  BRW  coordinate  system  [V)|. 

The  coordinate  transformation  from  [V],,  to  [V]^  is 
given  by 

tV3„  - h.(86)  [V],  (S-Ioia) 

where  R,(9j)  is  the  rotation  operator  which  rotates  the  [V], 
system  by  the  rotation  of  bisector  angle  about  Using 

the  explicit  expression  of  rotation  operator.  Eg.  5-lOla  can 


(S-lOlb) 


The  coordinate  transformation  from  [V],  to  [V]  is  given  by 

m,  - R.(»|)  [7],  (5-3a) 

where  is  the  rotation  operator  which  rotates  the  [V] 

system  by  the  rotation  of  the  turntable,  *j  about  the  x^. 

The  coordinate  transformation  from  [V]  to  [V),  is  given  by 

[V]p  - IV],  - [*,],  (5-2a) 

The  transformation  from  [V],,  to  [V],  is  thereforei 

[V]„  = 6,(0^)  R>(«,)  aV],  - X,)  (5-102) 

where  rotation  operators  and  translation  vectors  are  given 
by  Eqs.  S-lOlb,  S-3b  and  5-2b. 

A dose  distribution  for  multiple  arcs  is  obtained 
simply  by  adding  each  contribution  from  each  arc,  and  can  be 
represented  on  dose  grids  defined  (cartesian  or  oblique).  A 
dose  model  for  multiple  arc  is  expressed  in  the  form  given 
by 


(5-S 


Where  j is  an  arc  nunber  and  is  given  by  Eq.  9-96. 

Figures  S-22a  and  b give  montages  of  spatial  contours 
of  dose  distribution  Iron  the  dose  model  developed  on  major 
arc  planes  for  100’  and  180'  partial  rotation  arcs  with 
iaocenter  position  (0,0,0)  and  collimator  size  1 cm.  The 
isodose  shape  in  figures  5-22a  and  b is  in  good  agreement 
with  that  obtained  from  the  exact  dose  model  (Figures  S-22C 

Figures  5-23a,c  and  e give  montages  of  spatial  contours 
of  dose  distribution  from  a modified  cylindrical  dose  model 
on  three  orthogonal  planes  (axial,  sagittal  and  coronal)  for 
four  standard  arcs  (used  in  routine  stereotactic 
radiosurgery)  with  isocenter  position  (0,0,0)  and  collimator 
size  1 cm.  The  isodose  shape  in  figures  9-23a,  c and  e is 
in  good  agreement  with  that  obtained  from  Che  exact  dose 
model  (Figures  S-23b,  d and  f). 

Figure  5-24b  give  montages  of  spatial  contours  of  dose 
distribution  from  a modified  cylindrical  dcse  mcdel  on  the 
coronal  plane  for  eight  standard  arcs  with  two  Isocenters, 
(0,-0. 6,0)  and  (0,0. 6,0)  (four  arcs  about  each  isocenter) 
and  collimator  size  1 cm.  The  isodose  shape  in  figures 
S-24b  is  in  good  agreement  with  that  obtained  from  the  exact 
dose  model  (Figure  5-24a),  and  shows  more  accurate  dose 
distribution  than  those  from  dose  models  previously 


Figure  5-22.  laodose  distributions  calculated  by  loodified 
cylindrical  dose  model  (a,b)  and  exact  dose 
modal  (c,d)  tor  partial  rotational  arcs.  (a) 
and  (c)  represent  isodose  distributions  for 
100'  arc  (gantry  30*  to  130*).  (b)  and  (d) 

represent  isodose  distributions  for  180  are 
(gantry  -90*  to  90').  Isodose  lines  displayed 
are  from  90  to  10  in  decrements  10. 
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Figure  5-23.  Isodose  distributions  calculated  by  modified 
cylindrical  dose  model  {a,c,e}  and  exact  dose 
model  (bpd.f)  for  four  standard  arcs  (100’  arc 
with  turn-table  orientation  0",  4S',  -45",  and 
-90')  with  1 cm  collimator  sire.  The 
distributions  are  displayed  onto  the 
transverse  (a,b),  sagittal  (c,d)  and  coronal 
planes  (e,f).  Isodose  lines  displayed  are 
from  10  to  90  in  increments  10  and  normalised 
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developed  (l.e.  the  spherical  dose  model  in  Flg.5*24d  and 
the  cylindrical  dose  model  for  full  rotation  arcs  in 

As  a consequence,  the  programs  required  approximately 
20  seconds  to  calculate  dose  for  four  standard  aros  on  a 20 
X 20  matrix  (400  points)  using  the  modified  cylindrical  dose 
model  for  partial  rotational  arcs  and  an  additional  7 
minutes  to  perform  grid  generation  to  represent  isodose 
curves  with  a Zenith  386/16  computer  without  floating  point 
processor.  This  is  much  faster  (approximately  20  tines) 
than  the  exact  dose  model  to  calculate  dose  for  four 
standard  arcs.  Thus,  with  the  use  of  a simple  doss  model 
optimum  parameters  can  be  found  20  times  faster  than  with 
the  exact  dose  model,  which  makes  it  possible  to  perform 


real-time  optimisation  during  computer  dosimetry  treatnent 


CHAPTER  6 

BXPERIHEKTXL  QPTIMIZATIOK 


Since  LIHAC-besed  stereotactic  radiosurgery  uses 
multiple  noncoplanar  arcs,  three-dimensional  dose  evaluation 
and  many  beam  parameters,  a lengthy  computation  time  is 
required  to  optimise  even  the  simplest  case  by  a trial  and 
error  or  a computer  search  method.  The  basic  approach 
presented  in  this  paper  is  to  use  both  an  experimental 
optimization  and  an  analytic  optimization  to  solve  the 
problem  quickly  and  efficiently. 

An  important  problem  is  the  utilization  of  the  proper 
beam  parameters,  to  which  dose  shape  is  sensitive,  as  the 
optimum  variables  to  obtain  a favorabls  result  economically. 
In  comparing  optimum  variables,  multiple  isocenter  positions 
and  collimator  size  and  shape  are  most  useful  for  changing 
the  shape  and  margin  of  the  high  isodose  surface,  which 
makes  it  possible  to  shape  the  field  to  the  target  by 
varying  these  parameters.  Since  arc  variables  such  as 
length,  direction,  and  weighting  do  not  change  much  about  a 
high  isodose  shape  within  a restricted  range  of  arc 
variables  with  fixed  isocenter  positions  and  collimator 
sizss,  it  is  possible  to  change  the  low  isodose  shape 


encompassing  critical  organs  vtiila  maintaining  the  shape  of 
the  target  dose. 

Experimental  optimization  is  baaed  on  two  approaches. 
One  is  shaping  the  target  volumes  through  the  use  of 
multiple  isocenters  determined  from  dose  experience  and 
testing.  The  other  method  is  conformal  therapy  using  a 
beam's  eye  view  technique  and  field  shaping. 

Two  exact  dose  models  developed  for  circular  and 
rectangular  fields  in  chapter  5,  are  used  for  experimental 
or  visual  optimization.  The  circular  dose  model  is  used  for 
a multiple  isocenter  approach,  while  a rectangular  dose 
model  is  used  for  the  field  shaping  approach. 

Hultlple  Isoeentar  Approach 

In  this  section,  a brief  description  of  a multiple 
isocenter  approach  is  given,  and  a useful  rule  for  optimum 
isocenter  separation  and  collimator  size  is  developed  to 
shape  the  target  margin  uniformly  with  an  80%  isodose 
surface  for  an  elongated  target. 

The  method  is  based  on  the  relationship  between  an 
isodose  shape  and  optimum  variables.  The  current  technique 
utilizee  a single  isocenter  approach  with  multiple 
nonccplanar  arcs.  This  approach  results  in  spherical  dose 
distributions  in  the  target  and  dose  fall-offs  outsids  ths 
target,  which  depend  on  the  arrangement  of  the  arc  system. 
The  method  presented  here  is  to  use  multiple  isocenters  with 


standard  arcs  to  shape  target  volunes  through  the  use  of 
multiple  spherical  targets.  The  approximate  spherical  dose 
model  may  be  used  to  see  the  dose  shape  imaediately. 

The  test  is  based  on  the  use  of  two  or  three  isocenters 
in  parallel  for  l-D  shape  about  an  elongated  target.  Dose 
distributions  were  inspected  for  different  isocenter 
separations  and  collimator  sizes.  Four  standard  arcs  (three 
100*  and  one  ISO')  with  equal  arc  spacing  (45',370*  315', 
and  0*  turntable  angle)  were  used  for  each  isocanter.  After 
checicing  the  dose  shape,  including  field  uniformity  and  dose 
falloff,  the  useful  combinations  of  isocenter  separations 
and  collimator  sizes  were  considered. 

Figures  6-la  to  f represent  two  major  beam  profiles  for 
three  different  collimator  sizes  and  different  isocenter 
separations  with  two  isocenters:  horizontal  axes  along  two 
isocenters  on  the  coronal  plane;  and  vertical  axes 
perpendicular  to  horizontal  axes  at  the  mid-point  of  two 
isocenters  on  the  coronal  plane.  Figures  6-2a  to  d 
represent  isodose  distributions  on  the  coronal  plane  for 
different  isocenter  separations  with  1 cm  collimator  size. 
Figures  6-3a  to  d show  3-D  dose  profiles  on  the  sane  coronal 

A reasonable  optimization  criteria  for  the  target  dose 
is  to  fit  the  target  margin  within  the  sot  isodose  surface. 
When  we  use  multiple  isocenters  in  parallel,  the  isodose 
surface  can  be  elongated-shape.  The  elongated  target  shape 


Figure  6-1.  Dose  profiles  through  horiiontal  axis  (a,o,e) 
and  vertical  axis  (b,d,f)  for  three  different 
colliioatot  sizes  and  four  different  iaocenter 
separation  using  two  isocenters.  The  numbers 
shown  in  figure  represent  the  isocenter 
separation  in  cm.  The  relative  dose  is 
normalized  at  isocenter. 


Figure  e-2.  Isodose  distributions  on  the  coronal  plane  with 
1 cm  collimator  site  for  different  isocenter 
separation:  (a)  0 cm,  (b)  0.4  cn.  (c)  1 cm,  (d) 
1.2  cm.  The  isodose  lines  displayed  are  from 
90  to  10  in  10  decrements  and  normalized  by 


Figure  6-3.  3-D  dose  profiles  on  Che  coronal  plane  for 

different  Isooenter  separation  in  figure  6-Z. 


peanut 


can  be  approxiaated  by  a aliptical,  a cylindrical  or 
shape  with  the  longest  and  shortest  nargin.  Figure  6-4 
represents  longest  (L)  and  shortest  distances  (S)  of  the 
target,  which  can  be  fitted  to  the  80%  isodose  surface  for 
different  isocenter  separations  and  collimator  sizes. 
Although  other  isodose  lines  can  be  used  to  fit  the  target 
volume  according  to  the  shape  of  the  dose,  the  guideline  for 
80%  may  give  some  basic  idea  to  obtain  the  best  target 
corresponding  dose.  These  guidelines  can  be  used  to 
determine  the  isocenter  positions  and  collimator  sizes  for 
the  target  of  interest. 

Figures  6-5a  to  f represent  two  major  bean  profiles  for 
three  different  collimator  sizes  and  different  isocenter 
separations  with  three  isocenters  in  parallel!  horizontal 
axes  along  three  isoeenters  on  the  coronal  plane;  and 
vertical  axes  perpendicular  to  horizontal  axes  at 
mid-isocenter  of  three  isocenters  on  coronal  plane.  Figures 

5- 6a  to  i represent  isodose  distributions  on  axial,  sagittal 
and  coronal  planes  for  the  optimum  isocenter  separations, 
which  are  60%  of  collinator  sizes  {1,2, and  3 cm). 

Figure  S-7  represents  Che  guidelines  for  three 
isocenters.  For  a more  elongated  target,  ws  can  add  another 
isocenter  to  obtain  the  desired  height  and  diameter  (Fig. 

6- 8a) . It  is  possible  to  use  more  than  three  isocenters  in 
parallel,  however  there  is  no  increased  benefit  if  the 


isccentef  Seoeration  (cm) 


Guideline  for  two  isocenters.  The  solid  lines 
represent  the  longest  margin  (L),  and  dashed 
lines  represent  the  shortest  margin  (S)  of  the 
elongated  target,  which  can  be  fitted  to  the  ao 
i isodose  surface. 


Dose  profiles  through  horizontal  axis  {a,c,e) 
and  vertical  axis  (b,d,£)  for  different 
collimator  size  and  isocenter  separation  using 
three  isocenters  in  parallel. 
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Isodose  distribution  on  the  three  orthogonal 
planes  (axial,  sagittal,  coronal]  using  three 
isocenters  in  parallel  and  the  optimum 
isocenter  separations  to  give  uniform  target 
dose,  which  are  60  % of  collimator  size. 

(a,b,c)  for  1 cm,  (d,e,f)  for  2 cm,  and  (g,h,i) 
for  3 cm  collimator  size.  The  isodose  lines 
displayed  are  from  10  to  90  in  increment  10  and 


(continued) 


f 


Isocenter  Seperation  (cm) 


Guideline  for  Uhree  isocenters  in  parallel. 
The  solid  lines  represent  the  longest  nargin 
(L),  and  dashed  lines  represent  the  shortest 
margin  (S)  of  the  assumed  elongated  target, 
which  can  be  fitted  to  the  80  t isodose 


Figure  6-8. 


Isodose  distribution  using  more  than  two 
isocenters  in  parallel.  (a)  three  isocenters 
with  collimator  size  of  1 cm  and  isooenter 
separation  of  1.2  cm,  (b)  five  isocenters  with 
collimator  size  of  1 cm  and  isocenter 

displayed  are  from  10  to  90  in  increment  10  and 
normalized  by  maximum. 
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r in  3-0  space.  Figure  6-9a  and  f 
*tributions  on  coplanar  planes  through 
iree  isocenters,  which  forms  a triangular  geometry  and 
Ives  uniform  target  dose  with  two  possible  isocenter 
sparations.  This  isocenter  separation  configuration  will 
! a basic  unit  to  fit  the  2-0  shape  of  a target  by 
rranging  or  adding  isocenters.  Figures  6-loa  and  f 
ipresent  dose  distributions  on  the  plane  through  four 
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study  of  multiple  isocenters. 
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find  the  optimum  isocer 
trial  and  error  method 


il  dose  model  is  much  faster  than 
j model,  it  would  be  very  useful  1 
;er  position  and  collimator  i 
ir  computer-aided  design  opt: 


Isodose  distribution  on  the  three  orthogonal 
planes  (coronal,  sagittal,  axial)  using  three 
isocenters  in  triangle  geometry  on  the  coronal 
plane  with  1 cm  collimator  size  and  two 

and  1.2  cm  (d,e,f).  The  isodose  lines 
displayed  are  from  10  to  90  in  increment  10  and 
normalized  by  maximum. 


Figure 


6-10. 


Isodose  distribution  on  the  three  orthogonal 
planes  (axial,  sagittal,  coronal)  using  four 
isocenters  in  tetrahedra^geomet^  uito^l^CB 

separations;  1 cm  (a,b,c)  and  1.2  cm  (d,e,f). 
•as  isodose  lines  displayed 


c 


technique.  The  use  of  computer-aided  design  optimization 
with  a spherical  dose  model  is  discussed  in  chapter  7. 
Another  possible  study  is  to  utilize  conformal  therapy  using 
beam's  eye  view  and  field  shaping.  This  technique  is 

In  conclusion,  the  guidelines  in  Figs.  6-4  and  6-7  are 
useful  in  determining  the  optimum  isocenter  position  and 
collimator  size  for  elongated  target  shapes,  and  much  more 
benefit  is  obtained  from  the  multiple  isocenter  approach 
rather  chan  the  single  isocenter  approach  to  fit  the 
elongated  target  shape  (e.g.  within  an  flO%  isodose  surface, 
compare  Fig.  6-lla,b  with  c,d) . 


In  designing  radiation  treatment  plans,  the  objective 
is  to  conform  and  confine  the  high  dose  regions  to  the  shape 
of  the  target  as  closely  as  possible.  It  is  necessary  to 
shape  the  aperture  of  the  beam  and  compensate  for  surface 
irregularities,  internal  inhomogeneities,  and  target  shape 
to  achieve  desirable  dose  distributions.  Since  radiosurgery 

compensation  problem  is  not  important.  Our  concern  in  this 
section  is  field  shaping  using  the  beam's  eye  view 
projection  technique  and  adjusting  of  the  collimator 


(rotation  and 


of  using  irreguls 


shap^target^aength  - 

“Sir 
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shapes  such  as  multilaaf  collimators,  the  variation  of  a 
rectangular  field,  which  seems  to  be  more  practical  in 
radiosurgery,  is  considered  in  this  section. 


For  practical  3-D  treatment  planning,  one  of  the  most 
important  perspectives  to  be  used  is  the  Beam's  Eye  view 
(BEV) , since  it  displays  the  relationship  of  the  target 
volume  to  the  proposed  radiation  beans.  Beam's  eye  view  is 
a display  of  relevant  anatomical,  dosimetric  and  beam  data 
from  the  perspective  of  the  radiation  source,  looking  down 
the  radiation  beam,  Discussion  of  the  beam's  eye  view 
concept  and  it's  implementation  at  other  sites  can  be  found 
in  the  literature  (Fraass  et  al.,  I987j  Soitein  et  al., 

19821  Purdy  et  al.,  1987;  Chin  et  al.,  1985;  Mohan  et  al., 
1988).  The  patient  structures  which  have  been  defined  can 
be  made  into  wire  frame  (Houlard  t Dutrelx,  1984)  or  shaded 
3-D  views  (Bauer-Kirpes  et  al.,  1987).  In  this  worX,  BEV  is 
accomplished  by  displaying  a wire-frame  diagram  of  contours 
outlining  the  patient's  anatomic  structures  in  three 
dimensions  as  if  being  viewed  from  the  source  of  radiation 
along  the  central  axis  of  the  beam. 

In  order  to  display  the  defined  anatomy  in  beam's  eye 
view,  transformation  was  performed  from  the  anatomical 
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coordinates  were  transformed  into  screen  coordinates.  The 
first  Is  called  viewing  transformation,  and  second, 
perspective  transformation. 

The  viewing  coordinate  system  is  affixed  to  the 
collimator  of  the  treatment  machine  with  its  x -axis 
pointing  from  the  isocenter  to  the  source  and  and  a^-axes 
parallel  to  the  jaws  defining  the  field  size  as  discussed  in 
rectangular  field  in  chapter  5.  The  origin  of  this 
collimator  system  is  at  isocenter.  The  contour  points 
defining  the  structures  are  given  in  the  reference  BRW 
coordinate  system,  which  are  transformed  from  the  original 
coordinate  system  of  Che  CT  images.  These  points  undergo  a 
series  of  rotational  and  translational  transformations  to  be 
expressed  in  terms  of  the  collimator  system.  From  Egs.  S-70 
to  5-73,  [VJ|  can  be  transformed  to  the  system  by  the 
combined  transformation: 


[Vic  - M^i)  R,(«|)  I EV],  - [I,],  ) (5-70) 


The  collimator  coordinates  can  be  expressed  in  more  detail 
in  terms  of  BEIH  coordinates 

X.  - (X,  - X,)  cos  e,  + [(y,  - Yj)  cos  dj  + (2,  - Z,) 

sin  dj  ] sin  e,  <5-75) 

y,  - <-(*,  - Xj)  sin  e,  ♦ ( [y,  - Yj)  COB  d,  •*- 

(2,  - Zj)  sin  dj)  cos  e,]  cos  + <-(y,  - Yj) 
sin  dj  + (z,-  Zj)  cos  dj  ) sin  r,  (5-76) 

*c  “ - C-(x,  - X|)  sin  9,  -I-  I (y,  - Y.)  cos  dj  + 

(z,  - Zj)  sin  dj)  cos  e,]  sin  t,  + (-(y,  - Yj) 

After  the  viewing  transformation,  those  points  are  then 
projected  on  the  defined  imager  screen.  Suppose  a point  D 
is  chosen,  whose  collimator  coordinates  are  (d,0,0).  The 
plane  x^’sd  Is  the  screen  perpendicular  to  the  x^  axis.  The 
screen  coordinate  system  is  then  determined  by  prescribing 
its  origin  to  be  D and  its  y,  and  z,-axes  to  have  the  same 
directions  as  the  x^  and  y,,-axes,  respectively.  To  take 
into  account  the  beam  divergence,  the  objects  nearer  to  the 
source  than  screen  are  proportionally  magnified  and  those 
farther  than  the  screen  are  demagnified  (Aaimeraal,  1987). 

In  our  implementation  of  the  beam's  eye  view,  the  isocenter 


the  sane  way  we  can  de 


(SAD  - X,) 


The  screen  distance  d can  be  any  distance,  however  tt 
isocenter  can  be  assumed  to  be  at  the  center  of  the  scr 
in  practical  beam's  eye  view.  A note  advanced  display 
method  could  also  be  used  with  the  beam  visualization 
algorithm  (Reynolds  et  al.,  1988).  However,  the  use  of 
algorithm  is  more  complicated  and  is  not  necessary  Cor 
study  purpose.  Using  Eqs.  5-75  to  77,  and  Eqs.  6-3  and 
the  program  was  written  in  BASIC  language  to  display  th 


Figure  6-12,  The 


Optimization  Using  BEV  and  Field  Shaping 


The  basic  idea  of  BEV  is  to  determine  the  optimum 
directions  and  field  sizes  of  oblique  noncoplanar  beams  so 
as  to  maximize  the  coverage  of  the  target  region  while 
minimizing  the  inclusion  of  healthy  critical  organs  in  the 
high  dose  region.  The  best  arc  direction  in  radiosurgery 
can  be  determined  among  the  possible  directions  by  observing 
the  dose  distribution  through  the  entire  volume  in  the 
patient.  Another  way  to  determine  the  arc  position  is  to 

volume  matrix  (Kyrianthopoulos  et  al.)  through  the 
information  obtained  from  beam's  eye  view.  A dose  volume 

provide  a complete  summary  of  the  entire  1-D  dose  matrix 
relative  to  the  irradiated  and  unirradiated  volume  of  the 
patient.  This  is  especially  helpful  when  evaluating  how 
much  of  a normal  organ  is  receiving  an  amount  of  dose.  It 
is  invaluable  when  several  plans  have  to  be  evaluated  and 

to  the  smallest  size  that  covers  the  target  volume.  The 
outline  of  the  aperture  should  conform  to  the  shape  of  the 
target  visible  from  the  direction  of  incidence  and  provide 
adequate  margins. 

In  practical  cases,  it  is  convenient  to  automatically 
specify  the  shape  of  the  aperture  enclosing  the  target 


This  is  especially  inporte 


for  moving 


region. 

treatments  using  multi-leaf  collimators.  In  this  study,  a 
rectangular  collimator  is  used  to  fit  the  target  margin 

rectangular  collimator  is  easy,  less  time  consuming, 
provides  a simpler  dose  model,  and  quite  accurate  for  a 

Figures  s-13  and  6-14  show  the  projection  of  a 

rectangular  and  cylindrical  targets  as  seen  from  the  source. 
The  targets  have  a 3-D  description  including  7 slices  using 
a wire-frame  model.  The  center  of  a target  is  located  at 
isocenter  in  the  spherical  head  model.  The  dimension  of  the 
target  and  all  other  geometry  are  shown  in  Figure  6-15.  The 
objective  in  each  view  is  changed  as  the  collimator  rotates 
(at  a ten  degree  angular  step)  with  four  different  table 
orientations  (O,  45,  -45,  -90).  The  fields  in  Figures  6-13 
and  6-14  are  selected  to  encompass  the  target  plus  some 
margin,  and  are  shown  as  dashed  lines.  Dose  distributions 
are  then  calculated  using  the  rectangular  dose  model 
developed.  Figures  6-16  to  6-17  show  three  orthogonal  dose 
distributions  for  cylindrical  targets  using  one  arc  (table 
angle  0 or  -90)  or  two  arcs  (table  angle  0 and  -90).  Dose 
distributions  in  figures  6-17a,  b,  and  c show  results 
comparable  with  that  obtained  from  the  multiple  isocenter 
approach  (d,e,f  in  Fig.  6-17).  Figure  6-18 


shows  thre 


Beam's  eye  view  of  a rectangular  target  in 
Fig.  6-15.  The  target  has  a 3-0  description 
including  7 slices  using  a wire-frame  model, 
(a)  table  orientation  O’,  (b)  -90’,  (c)  4S', 
and  (d)  -46’.  The  two  cross  lines  represent 
two  central  axis  of  O'  collimator.  The  dashed 
line  represents  rectangular  field  outline 
which  can  shape  target  projection  with  the 
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Figure  6-13  (cc 


Seam's  eye  view  of  a cylindrical  target  in 
Fig.  6-15.  The  target  has  a 3-D  descriptio 
including  7 slices  using  a wire-frame  model 
(a)  table  orientation  O’,  (b)  -90',  (c)  45’ 
and  (d)  -45’. 


'a.  e Gfi:  38 

Ifi:  -58  G?:  -138 

--4- 

1 

!-;•  i'  G8-  55 

Ifl:  -30  G^  -105 

Tft:  0 GR:  80 

— W — 

T8:  -38  Gll:  -88 

T 

Ii;-:  8 G8:  185 

Tfl:  -80  Gft:  -55 

T 

8 Gfl--  130 

Tft:  -58  Gfl:  -38 

A 


19:  45  G'S:  38 

I»:  -IS  fi;  -136 

iil:  45  G9:  55 

19: -45  G«  -5?5 

Ifi:  45  Sfl:  S8 

ih:  -43  Wi:  -HE 

Tii:  « Gft:  IK  T9:  -45'  Gfi:  -55 


-%r- 

45  ffl;  13E 

Ifi:  -45  G4:  -3f 

-4»^ 

r 

c 

D 

(continued) 


X 


B 


Figure  6-l«. 


Isodose  curves  on  the  three  orthogonal  planes 
(ajcial,  sagittal,  coronal)  through  the  target 
point  due  to  the  field  shaping  technique  using 

table  orientation  O'  with  gantry  angle  (-80* 
to  80')  or  (b)  -90'  with  gantry  angle  (30*  to 
130').  The  isodose  lines  displayed  are  from 
10  to  90  in  increment  10  and  normalized  by 


Isodose  curves  on  the  three  orthogonal  planes 
(axial,  sagittal,  coronal)  through  the  target 
point  due  to  (a)  the  field  shaping  technique 
using  beam’s  eye  view  from  two  arcs  at 
table  orientation  O’  with  gantry  angle  (-80’ 

130’),  and  (b)  due  to  two  multiple  isocenter 
technique  using  standard  eight  arcs.  Both  80 

the  cylindrical  target.  The  isodose  lines 
displayed  are  from  10  to  90  in  increment  10 
and  normalized  by  maximum 


orthogonal  dose  distributions  for  a rectangular  target  using 
angle  0,  45,  •45,  and  -90).  The  collinator  was  rotated  to 
achieve  proper  field  batching.  Figure  6*19  shows  a dose 
volume  histogram  for  a rectangular  target  in  two  cases  (4 
arc  field  shaping  and  2 multiple  isocenterj . 


Isodose  curves  an  the  three  orthogonal  planes 
(axial,  sagittal,  coronal)  through  the  center 
of  rectangula  target  due  to  the  field 
shaping  technique  using  beam's  eye  view  from 
four  arcs  at  table  orientation  0*  with  gantry 
angle  (-80'  to  80'),  -90',  45’  and  -45"  with 
gantry  angle  (30’  to  130') . The  isodose  lines 
displayed  are  from  XO  to  90  in  increment  10 
and  normalized  by  maximum. 
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6-19.  Dose  VoliBoe  Histogram  for  a 


rectangular  target 


CHAPTER  7 

ANALYTICAL  OPTIMIZATION 


coaputar-Alded  Paaign  ODtlnlzatlon 

A computer-aided  design  (CAD)  optimization  is  any  use 
of  an  algorlthnlc  approach  to  optimizing  a dose  distribution 
using  a quantitative  evaluation  of  the  optimality  of  a 
treatment  plan,  in  this  category  of  optimization  methods, 
one  expresses  the  optimization  criteria  as  some  mathematical 

evaluates  the  relationship  of  the  calculated  dose 


the 


since  it  is  not  efficient  to  quantify  all  the 
opclnization  vatiables  in  CAD  optimization,  the  better 
approach  is  to  utilize  the  most  sensitive  beam  parameters  as 


From  the  study  of  the  effects  of  optimum  variables, 
isocenter  positions  and  collimator  sizes  are  useful 


varying 


sighting 


target  dose.  A cylindrical  dose  aodel  with  one  full 
rotational  or  partial  rotational  arc  is  used  with  a 
computer-aided  design  optimization  technique.  The  third 

exact  dose  model.  The  isocenter  position  and  collimator 

selected  by  the  treatment  planner,  leaving  the  set  of  arc 
parameters  as  the  variables  which  are  left  for  the  CAD 
optimization.  The  objective  function  is  then  written  as 
some  function  of  the  arc  variables,  and  optimized  through 
CAD  optimization  technique. 

In  the  present  work,  the  objective  function  is  selected 
to  minimize  the  weighted  dose  to  the  critical  organs  defined 
at  different  positions,  or  to  maximize  the  dose  gradient 
between  target  boundary  and  normal  structure  around  target. 
He  can  now  write  our  nonlinear  constrained  optimization 
problem  mathematically  as  follows: 

Minimize:  r Dj(x)  Objective  function  (7-ia) 

Maximize:  £ (D,  (s) -D„(x)  )■  (7-lb) 

Subject  to: 

d'  s D,(»)  S D"  constraints  (7-2) 


snstraints 


'-3) 


Design  vector 

Dose  to  the  critical  organ  a 


Lower  and  upper  limits  c 


t position  t 

a position  n 

e design 


A more  general  category  of  algorithms,  referred  to  as 
nonlinear  programming  (NLP) , is  needed  to  solve  this  general 
optimization  problem  (Jacoby  et  al.,  1972;  Ben-Israel  et 
al..  19S1;  Shoup  i Mistree,  1987]. 

I^early  all  of  the  optimization  techniques  that  are 
currently  in  use  can  be  categorized  as  methods  for 
unconstrained  minimization  or  maximization.  This  means  that 
they  find  an  optimum  without  regard  to  any  constraints. 

Since  our  problems  have  a number  of  equality  and  inequality 
constraints,  the  use  of  an  unconstrained  optimization 
algorithm  by  itself  will  usually  not  provide  the  best 
approach  to  finding  a usable  solution.  The  penalty  function 
approach  is  one  powerful  way  to  adapt  constrained 
optimization  problems  for  solution  by  unconstrained  methods 
(Fiacco  & McCormick,  1968] . 


Ity  function 
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ipolation  along  a lit 


3en  the 


decreased  below  a predetermined  value  with  no  substantial 
change  in  the  objective  function,  the  procedure  algorithm 

elegance,  it  is  a highly  efficient  optiniaation  algorithn. 
The  relative  simplicity  cf  the  method  has  made  it  one  of  the 
most  popular  for  computer  implementation. 

Powell  Method 

to  it's  original  development,  is  one  of  the  most  efficient 
and  reliable,  and  certainly  one  of  the  most  popular  of  the 

of  conjugated  directions,  where  directions  S'  and  s'  are 
conjugate  if 


in  n orthogonal  directions  (i.e.  n design  variables),  s', 
i-l,n  being  the  coordinate  directions,  where  each  search 


(7-S) 


where  q is  the  iteration  number,  x is  the  vector  of  design 
variables,  B°  is  the  search  direction,  and  is  a scalar 
multiplier  determining  the  amount  of  change  in  x for  q 
iteration.  Having  completed  the  n unidirectional  searches, 
a new  search  direction  is  created  by  connecting  the  first 
and  last  design  points.  This  becomes  the  n+1  search 
direction.  It  is  convenient  to  store  the  search  information 
in  a Hessian  array  matrix  R.  He  begin  with  the  identity 


(7-6) 


The  columns  of  R are  the  unidirectional  search  vectors  s', 
iei,n.  After  finding  the  minimum  in  each  direction  (i.e. 
one  dimensional  search  using  golden  section  or  quadratic 
interpolation),  we  replace  s'  in  matrix  H by  a, s'  to  give, 
after  n iterations: 


conjugate  direction: 


>rch 


dire 


column  n.  This  provides  a new  B natrix  containing  n search 


algorithm  for  Powell's  nethod  is  given  in  Figure  7-2. 


Penalty  Function  Method 

A conmon  approach  to  solving  the  constrained 
minimization  problem  is  based  on  the  penalty  function 
approach  with  the  unconstrained  minimization  nethod.  The 
main  idea  of  the  penalty  function  approach  is  to  minimize 
the  objective  function  as  an  unconstrained  function  but  to 

critical  organ  is  outside  the  constrained  arc  range,  a 
certain  penalty  is  assessed  to  the  objective  function.  The 

create  a pseudo-objective  function  by  summing  the  penalty 
function  and  the  previous  objective  function: 


where  F(»)  is  the  original  objective  function,  which  is 
represented  by  Eg.  7-1.  P(s)  is  an  imposed  penalty 

ch  depends  on  the  penalty  function 


function,  the  form  of  whic 


Figure  7*2.  The  Potfell  conjugated  direction  search 
algorithm. 


sing  ample 


ilty  functic 


penalty  function;  or  an  extended  interior  penalty  function. 
The  penalty  function  is  zero  for  all  design  points 
satisfying  the  constraints.  Clearly,  if  all  constraints  are 
satisfied,  the  mininization  of  »(*,r  ) is  equivalent  to  the 
mininization  of  F(x).  The  scalar  r is  a multiplier  which 
detenaines  the  magnitude  of  the  penalty,  and  is  held 
constant  for  a complete  unconstrained  minimization.  The 
subscript  p is  the  unconstrained  mininization  number.  The 
vaotor  X represents  the  optimum  variables. 

Rule-Based  Search  Method 

"conjugated  direction  search"  methods  are  well  suited  to 
finding  optimum  arc  positions  using  the  cylindrical  dose 

finding  optimum  isocenter  positions  and  collimator  sizes 
with  the  spherical  dose  model  because  they  find  local 
optima  instead  of  global  optima.  A simplistic  approach  is 
to  vary  the  relevant  parameters  in  suitable  increments,  and 

objective  function.  Since  the  number  of  variables  is  too 
large  to  compute  the  radiation  isodoses  for  all  passible 
configurations,  a total  search  method  is  not  desirable, 
consider,  for  exam 


iiple,  the  very  simple 


vary  as  follows  for  each  arc:  (1)  table  angle  0 • 180*  in 
10*  increnents  (IS  steps);  (2)  collimator  sizes  1.0  - 3.0  cm 
in  0.2  cm  increments  (11  steps);  (3)  isocenter  positions  in 
3D  space  5 x 5 x 5 mm  cubes  in  1 mm  increments  (125  steps) ; 
(4)  relative  arc  weighting  from  0.2  to  l.o  in  0.2  increments 
(5  steps).  The  number  of  choices  for  each  arc  is  than 
123,750  (i.e.  13  X 11  X 125  X 5).  The  total  number  of  plans 
to  be  calculated  is  then  123, 750^  or  2.345  X 10^.  The 
calculation  of  a single  plan  with  four  arcs  take  2 minutes 
to  compute  a 400  point  dose  matrix  using  a 386-16  computer, 
even  with  a math  coprocessor.  Then,  8.82  x 10^^  years  would 
be  required  to  run  through  the  various  options.  Even  the 
fast  approximate  dose  model  can  only  speed  computations  by  a 
factor  of  100.  Even  though  the  "pattern  search"  and 
"conjugated  direction  search"  methods  are  more  sophiscated 
than  the  total  search  method,  these  two  search  methods 
sometimes  have  problems  in  finding  global  minima  or  maxima. 

h rule-based  search  algorithm  was  developed  to  find  the 
optimum  isocenter  positions  and  collimator  sizes  to  maximize 
the  dose  gradient  between  the  target  boundary  and  normal 
structures  while  overcoming  the  local  minimum  problem,  one 
or  two  isocenters  are  considered  in  the  present  study,  since 
the  dose  gradient  using  more  than  two  isocenters  is  not 
sharper  than  that  using  a single  isocenter  in  fitting  the 
dose  distribution  to  a complicated  target  shape. 


The  use  of  this  search  nethod  consists  of  tvo  aajor 
steps.  First  is  the  determination  of  reasonahle  initial 
variables  from  the  imaged  target  shape,  and  second  is  to 
follow  a rule-based  search  procedure  according  to  the 
relationship  between  isocenter  separation  and  collimator 
sizes,  while  satisfying  constraint  conditions.  Prom  the 
patient  data  obtained  from  a series  of  CT  or  KRI  slices,  the 
target  can  be  assumed  to  be  a spherical  shape  or  can  be 
separated  into  tvo  spherical  targets  using  the  same  or 
different  collimator  sizes.  The  one  isocenter  case  is 
simpler  than  the  two  isocenter  approach,  so  we  will  discuss 
two  isocenters.  The  first  step  is  to  determine 
approximately  the  initial  beam  parameters  (isocenter 
positions  and  collimator  sizes)  from  the  3-D  target  shape 
which  is  generated  from  a series  of  target  contours.  The 
collimator  sizes  are  initially  chosen  to  be  smaller  than  the 
estimated  target  size,  so  that  the  target  volume  will  be 
included  in  the  planner  selected  isodose  surface  only  by 
increasing  the  collimator  size.  After  initial  values  are 
given,  the  second  step  is  to  determine  a search  vector  to 
efficiently  find  the  optimum  Isocenter  positions  and 
collimator  size  to  maximize  the  dose  gradient  around  the 
target,  instead  of  considering  an  arbitrary  search 
direction  (isocenter  positions  in  arbitrary  3-D  target 
area),  tha  direction  of  a vector  through  the  two  initial 
isocenter  positions  or  two 


arbitrary  positions  close 


initial  isocenter  positions  could  be  useful  assuoing 


initial  values  are  reasonable.  The  third  step  is  to  Chech 
the  objective  function  at  every  movement  of  the  isocenter 
position  in  a 1-D  search  direction  determined  in  step  2 as 
well  as  the  variation  of  collimator  size.  From  the 
relationship  between  optimum  isocenter  separation  and 
collimator  size  for  uniform  target  dose,  isocenter 
separation  can  be  adjusted  by  choosing  two  collimator  sizes 
(i.e.  optimum  isocenter  separation  ~ 80  to  120  percent  of 
Che  average  of  the  two  collimator  sizes),  and  two  temporary 
optimum  isocenter  positions  can  be  determined  by  moving  the 
initial  isocenter  positions  according  to  the  propotional 
ratio  of  collimator  sizes.  Therefore,  instead  of  checking 
every  point  on  the  search  vector  determined,  a couple  of 
points  checked  near  the  ekperimental  optimum  isocenters  is 
sufficient.  The  fourth  step  is  to  check  whether  constraint 
conditions  for  the  target  are  satisfied.  All  selected 
constraint  points  along  the  target  outline  are  checked  to 
ensure  that  the  planner  selected  isodose  surface  covers  all 
the  constraint  points.  One  of  two  collimator  sizes  is 
increased  according  to  the  comparision  of  isocenter-to- 
target  constraint  distance  normalized  by  collimator  size 
until  all  the  target  constraints  are  satisfied.  Than,  the 
isocenter  positions  are  readjusted  by  the  change  of 
collimator  sizes.  The  value  determined  here  ie  optimum  for 
the  first  iteration.  The  fifth  step  is  to  choose  another 


first  optimum  isocenter 


isocenter  positions  near  the 
position  determined  on  the  first  iteration  to  check  whether 
a better  optimum  can  be  obtained,  while  fixing  the  second 
optimum  isocenter  position.  The  procedure  is  repeated 
following  steps  one  through  four.  The  final  step  is  again 
to  fix  the  first  optim^Im  isocenter  position  while  varying 
the  second  iaocenter  position  near  the  initial  second 
isocenter  position.  The  change  between  first  and  second 
isocenter  is  repeated  until  the  planner  is  satisfied. 
However,  after  better  optimum  values  are  obtained  after 
several  cycles,  it  might  be  better  to  check  iscdose 
distributions  visually  using  the  exact  dose  model,  and  then 
readjust  slightly.  The  steps  for  this  method  are  given  in 

Aopllcatlon  and  Results 

He  discussed  two  main  approaches  and  detailed  methods 
to  optimize  dose  distribution  for  stereotactic  radiosurgery. 
The  implementation  of  computer-aided  design  optimization 
with  experimental  optimization  and  it*s  application  is 
discussed  to  illustrate  the  method  and  the  results. 

The  specific  problem  is  to  minimize  the  dose  to  the 
critical  organs,  which  are  represented  by  a few  points. 


while  giving  the  desired  dose  to  the  target  area.  The 
target  is  assumed  to  be  an  elongated  cylindrical  shape 
(height-1  cm,  diameter-0 . 5 cm),  and  is  located  at  the  center 
of  the  spherical  head  model  (20  cm  diameter).  The  critical 
organs  are  assumed  to  be  located  at  three  different 
positions  for  two  different  cases  (Figs.  7-4a  and  b] : case 
1-1,  e,(0,0,2),  Cj(0,2,0),  and  Cj(0,2,2):  case  1-2, 

C,(0,0,2),  02(0,2,0],  and  0,(0,2,0.828).  We  optimize  this 
problem  using  both  experimental  and  analytical  approaches. 

Fxteriinental  eotimization 

The  isocenter  positions  and  collimator  sizes  are 
predetermined  experimentally.  The  test  is  based  on  0-D 
(sphere) , 1-D  (elongated  shape) , 2-D  (plate  shape) , and  3-0 
(arbitrary  shape)  target  shapes.  Usually,  four  arcs  with 
egual  arc  spacing  (45*)  are  used  for  each  isocenter,  cases 
1-1  and  1-2  simply  deal  with  the  1-0  target  shape  problem. 
From  guidelines  shown  in  Pig.  6-4  developed  for  sphere  and 
elongated  target  shapes,  the  two  isocenter  positions 
(0, -0.5,0)  and  (0,0. 5,0)  with  collimator  size  1 cm  are 
chosen  for  case  1-1  and  2.  Now  the  80%  isodose  surface 
sufficiently  covers  the  target  shape  for  both  case  1-1  and 

Obiective  function  to  be  minimized 

After  the  isocenter  positions  and  collimator  sizes  are 
determined  experimentally  with  the  standard  arc  system,  the 


best  arrangafflsnt  of  arc  orientation  and  weighting  to 
sininize  the  doae  to  the  critical  organs  is  searched  for  by 
the  conputer-aided  design  optimization  technigue. 

In  the  present  worh,  the  objective  function  is  selected 
to  minimize  the  sum  of  the  weighted  dose  to  the  critical 
organs  defined  at  designated  positions.  Our  aim  is  now  to 
find  optimum  parameters  (arc  orientation  and  weighting)  such 
that  the  objective  function 


is  minimized  subject  to  the  constraint  conditions  described 
in  Eqs.  7-2  and  7-3,  where  c's  represent  positions  of 
critical  organs  designated  for  cases  1-1  and  1-2.  is 
relative  dose  to  the  critical  organ  located  at  position  c, 

given  by  Eg.  S-91.  All  the  arc  weights  are  assumed  to  be 
one.  The  critical  organ  weights  are  assumed  to  be  one 
in  case  1-1  , however,  much  more  weight  is  given  to  the 
position  Cj(0, 2 , 0. 228]  in  case  1-2. 


Constraints  chosen 

In  the  course  of  minimizing  the  sum  of  the  weighted 
doses,  the  constraints  (or  behavior  constraints)  and  side 
constraints  must  be  specified  and  quantified.  The  major 
constraint  condition  for  uniform  target  doss  is  to  assure 
the  delivery  to  the  target  volume  of  a relative  dose  larger 


certain  level  (e.g.  30%  or  70%),  which  is  expressed  by 
7-2.  The  side  constraints  require  a reasonable  range  of 
iables  due  to  physical  limitations  (e.g.  upper  and/or 
er  limits),  which  are  expressed  by  Eq.  (7-3). 

when  we  adjust  only  the  shift  angle  S of  the  total  arc 
entation  with  equal  arc  spacing  (4S'),  which  would  not 
nge  the  target  dose  shape,  the  oonstraint  conditions  are 


pllfled  to: 

0 S 1,  S T/%  (7-10) 

«,  - <,  (7-11) 

+ r/4  (7-12) 

Sj  = 4,  * T/2  (7-13) 

<(.  = «,  + 3T/i  (7-14) 

0 < 4j  S »/4  (7-15) 

#5  - (7-16) 


(7-17) 


where  end  6^  are  shift  angles  of  two  isocenters  froa  the 
standard  arc  position  and  0^  through  0^  are  table  angle 
orientations  for  eight  individual  arcs,  which  can  be 
detemined  easily  by  the  equality  constraint  relationship 
between  S and  0 (or  we  may  search  for  0's  directly  using 
different  constraint  conditions  and  a reasonable  arc  range] 

Penalty  function  method  for  constrained  nininization 

The  two  main  constraint  equations  (Eg,  ?-l0  and  7-lS) 
can  now  be  stated  in  the  standard  form  to  be  implemented  in 
the  penalty  function  approach: 


g,(a)  = - 

«,  < 0 (7-20) 

, - m/4  S 0 (7-21) 

gj(«)  - - 

< 0 (7-22) 

94(»)  = f 

j - m/4  < 0 (7-23) 

The  most  come 

in  penalty  function  used  in  Che  interior  method 

P(r)  = r 

(7-24) 

g,(«) 

where  P(z)  is  an  inposed  penalty  function.  The  vector  z 
represents  the  arc  variables  6^  and  6^  in  this  problem, 
other  arc  variables,  p,  through  Oj,  are  determined  easily 
iron  if  and  ig.  Our  pseudo-objective  function  is  now 
expressed  explicitly  by 

*(z,rj)  - F(z)  + r^  P(z)  (7-as} 

where  F[x)  is  the  original  objective  function,  which  is 
represented  by  Eq.  7-la.  He  optimize  this  problem  using  an 
unconstrained  minimization  algorithm  and  penalty  function 
method. 

Cases  1-1  and  1-2  were  solved  using  two  design 
variables,  and  ig,  by  the  interior  penalty  function.  For 
the  initial  values  of  arc  position,  a 2*  shift  angle  for  5f 
and  Sg,  which  is  close  to  standard  arc  position  (O'J,  was 
used.  In  each  case,  the  Hooke  and  Jeeve  method  and  the 
Powell  method  were  used  for  the  unconstrained  minimization 
subprogram.  The  initial  design  vector  and  optimum 

design  for  the  two  cases  are  given  in  table  7-1.  The  values 
obtained  from  the  computer  algorithm  seem  to  be  quite 
reasonable.  Initially  the  arc  positions  were  assigned  to 
hit  the  critical  organs.  The  final  arc  positions  were 
rearranged  not  to  hit  the  critical  organs  while  minimizing 
the  weighted  dose  to  the  critical  organa  in  both  case  1-1 
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n while  sacrificing  the  less  weighted 


critical  organ 


As  a consequence,  the  prograns  required  approximate!/  1 
minute  to  find  the  optimum  parameters  in  this  problem  using 
Che  Powell  method  with  the  full  rotational  cylindrical  dose 
model,  and  an  additional  2 min  for  the  Hooke  and  Jeeves 

processor.  This  speed  can  be  increased  with  a more  advanced 
computer  and  floating  point  processor,  which  makes  it 
possible  to  do  real-time  optimiaation  during  computer 
dosimetry  treatment  planning. 


iihlle  keeping  equal  arc 


spacing  and  may  give  a more  desirable  dose  distribution 
under  certain  conditions.  In  the  follow 


iiing  example 


consider  this  situation  and  discuss  specific  cases  to 
inpienent  the  CAD  optimization  technique. 

Problem  chosen 

Case  2-1:  The  target  is  assumed  to  be  a spherical 

shape  (tadiuB«0.8  cm)  located  at  the  center  of  the  spherical 
head  (diameter  - 20  cm).  The  critical  organs  are  located  at 
two  different  positions,  C, (0, -0 . 2 , 1 . 0)  and  Cj(0 , 0. 4 , -0 .8) 
in  the  reference  BRW  coordinate  system  (Pig.  7-5a). 

Case  2-2:  The  target  is  assumed  to  be  an  elongated 

cylindrical  shape  (height-2.0  cm,  diameter-0.4  cm)  located 
at  the  center  of  the  spherical  head  (dianeter-20  cm) . The 
critical  organs  are  located  at  two  different  positions, 
Cj(0,-o.6,l.0)  and  Cj(0,o.6,-1)  in  the  reference  BRW 
coordinate  system  (Fig.  7-5b). 

We  optimize  case  2-1  and  case  2-2  using  both  an 
experimental  approach  and  an  analytical  approach.  The 
isocenter  positions  and  collimator  sizes  are  predetermined 
to  shape  to  the  target  experimentally,  and  are  followed  by 
an  analytic  optimization  technique  to  search  for  the  best 
arc  parameters  (arc  spacing  or  arc  weighting)  to  minimize 
dose  at  the  critical  points  chosen. 

Experimental  optimization 

The  isocenter  positions  and  collimator  sizes  are 
predetermined  experimentally,  case  2-1  and  2-2  simply  dealt 
with  the  0-D  and  1-D  target  shape  problems.  From  guidelines 


VERT 


with  the  0-D  and  1-0  target  shape  problems.  Proa  guidelines 
shown  in  Pig.  6-4,  developed  for  spherical  and  elongated 
target  shapes,  the  single  isocenter  position  (0,0,0)  with 
collimator  size  1 on  is  chosen  for  case  2-i.  Two  isocenter 
positions  (0,-0. 6,0)  and  (0,0. 6.0)  with  collinator  site  1 cn 
are  chosen  for  case  2-2.  How  the  ao%  isodose  surface 
sufficiently  fits  the  target  shape  for  both  case  2-1  and 


Objective  function  to  be  minialyeH 

Our  aim  is  now  to  find  optimun  paraaeters  (arc 
orientation  and  weighting)  such  that  the  objective  function: 

* ° (7-la) 

is  minimised  subject  to  the  constraint  conditions  described 
in  Eqs.  7-2  and  7-3,  where  c's  represent  positions  of 
critical  organs  designated  for  case  3-1  and  case  2-2.  D is 
relative  dose  to  the  critical  organ  located  at  position  e, 
and  is  represented  by  the  approximate  dose  model  given  by 
Eg.  S-91.  All  the  arc  weights  are  assumed  to  be  one,  and 
the  critical  organ  weights  W,,  are  assiu&ed  to  be  one  in  both 


Constraints  chosen 

Instead  of  using  Eg.  7-2  as  constraints  to  maintain 
target  dose  shape,  it  may  be  simpler  and  more  desirable  to 


create  a few  constraint  conditions  that  can  search  for  the 
optiaal  result  by  varying  parameters  in  a restricted  range 
about  the  starting  point  values.  Initially,  four  arcs  with 
equal  spacing  (45'1  are  used  for  the  standard  plan.  From 
experience,  four  arcs  with  a small  variation  from  equal  arc 
spacing  (up  to  30')  does  not  significantly  change  Che  high 
isodose  shape.  Now  the  constraint  conditions  given  by  Eg. 
7-2  are  simplified  for  case  2-1  by 

-v/2  S 0,  S 0 

0j  > 0,  + r/6 

0j  i 0j  * »/6 

0j  * fr/6  < 0^  S T/3 
with  additional  conditions  for  case  2-2: 

-»/2  S 0,  S 0 (7-30) 

0j  2 0,  ♦ »/6  (7-31) 


(7-26) 

(7-27) 

(7-26) 

(7-29) 


0,  + <r/6  < 0j  £ r/3  (7-33) 

where  0,  to  0g  are  turntable  angle  orientations  for  eight 
arcs.  These  simplified  constraints  can  guic)lly  solve  the 
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treatment  planner  to  explore  unusual 
might  yield  better  results  than  more 
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by 


-ir/2 


T/2  . i-l,2...8 


Since  all  the  constraint  conditions  given  by  Eqs.  7-26  to 
7-33  satisfy  side  constraints  given  by  Eg.  7-34,  side 


Penalty  function  method  for  constrained  ninimlaatlBn 
shown  in  Eqs.  7-26  to  7-33  can  now  be  stated  in  the  standard 


g2(X)  = -a,  - r/2  S 0 (7-36) 
g,;*)  - a,  - Oj  t »/«  S 0 (7-37) 
g.(x)  = Oj  - Oj  + r/6  S 0 (7-38) 
g5(«)  = Oj  - 04  + »/6  S 0 (7-39) 
9j(»)  - 04  - r/3  S 0 (7-40) 


additional  conditions  for 


g,(a)  - *5  S 0 (7-41) 
g,(X)  "•  -«5  - r/2  s 0 (7-42) 
g«(*)  • «5  - *j  + ir/6  S 0 (7-43) 
9io<«)  = «6  - 07  * S 0 (7-44) 
g„(X)  = d,  - 7-  r/6  S 0 (7-4S) 
0i2(»)  = 0j  - T/3  < 0 (7-4S) 


The  penalty  function  used  in  the  interior  method  is: 


where  J * 6 for  case  2-1,  and  J * 12  for  case  2-2.  The 
vector  X represents  the  arc  variables  to  d&  in  case  1 and 
d|  to  dg  in  case  2.  Our  pseudo-objective  function  is  now 
expressed  explicitly  by  Eqs.  7-25,  7-la  and  7-24,  and  can  be 
minimized  through  the  unconstrained  optimization  algorithm 
programed. 

Feasible  starting  values  of  design  vector 

Because  of  the  nature  of  the  Inverse  function  from  the 
interior  penalty  function  method,  the  pseudo  objective 
function  »{x,r„)  will  approach  infinity  on  the  constraint 
boundaries.  The  solution  method  using  interior  penalty 


functions  proceeds  as  follows.  First,  a feasible  design 
point  in  the  allowable  space  is  chosen  inside  the 
constrained  set.  once  the  starting  point  is  inside  the 
constrained  area,  the  optimum  found  will  also  be  within  this 
area,  once  the  minimus  is  achieved,  it's  location  becomes 
the  starting  point  for  a new  pseudo-objective  function  in 
which  the  value  of  tp  is  reduced.  The  optimization  process 
is  repeated  for  a sequence  of  decreasing  tp  values. 

In  addition  to  Che  requirement  of  feasible  starting 
values,  good  estimated  starting  values  within  the  feasible 
range  saves  iteration  time  to  reach  the  optimum.  In  the 
present  problem  we  nay  choose  standard  arc  positions  as 
starting  values  of  arc  variables.  However,  when  Che 
standard  arc  hits  the  critical  organ,  it  is  desirable  to  use 
as  starting  values  Che  arc  positions  which  do  not  hit  the 
critical  organ.  This  speeds  reaching  for  the  optimum. 

Suppose  we  Xnow  the  coordinates  of  four  corners  of  the 
rectangle  which  sufficiently  covers  the  critical  organ  on  a 
series  of  CT  slices,  and  the  coordinate  of  the  isocenter 
position  in  the  BRH  coordinate  system.  The  purpose  is  to 
find  the  critical  arc  ranges  which  hit  the  critical  organ 
defined  by  a series  of  rectangles.  Using  the  geometricai 
relationship  between  spherical  and  cartesian  BRH  coordinate 
systems: 


(7-48) 


{7-49) 


R sin  s sin  « = - 


(7-SO) 


(7-51) 


(7-S3) 


A series  of  sets,  («S,e),  can  be  determined  for  the 
coordinates  of  the  rectangular  comers  from  Egs.  7-S2  and 
7-53.  Then,  the  critical  arc  range  (lower  and  upper  limits 
of  turntable  angle  o and  gantry  angle  9)  in  figure  7-6  are 
determined  from  this  information.  The  critical  arc  range 
(9'-9“  and  determined  can  be  constraint  conditions  or 

may  be  useful  to  determine  the  starting  arc  positions  in 
computer-aided  design  optimization. 

Case  2-1  and  case  2-2  were  solved  using  four  and  eight 
design  variables  by  the  interior  penalty  function  method. 
For  the  starting  values  of  arc  position,  turntable  angles 


Z, 


Figure  7-6.  Geooetrical  relationship  between  critical  arc 
range  and  isocentric  SRH  coordinate  eysten. 
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Figure  7-7a  gives  a nenCage  of  spatial  contours  of  dose 
distribution  on  the  coronal  plane  from  a modified 
cylindrical  dose  model  and  the  four  optimum  arc  positions 
searched  by  the  Powell  method  for  case  2-1.  The  isodose 
shape  in  figure  7-7a  shows  a comparable  result  with  that 
obtained  from  the  exact  dose  model  near  low  isodose  line 
(Figure  7-7b) . The  dose  distribution  obtained  with  final 
search  values  in  Fig.  7-7h  gives  improved  distribution 
compared  to  the  dose  distribution  obtained  with  the  starting 
design  in  Fig.  6-4a,  while  Keeping  the  target  margin  and 
minimizing  dose  near  the  critical  organs.  The  search 
doesn't  seen  to  be  complete  because  of  a rough  termination 
criteria  or  local  minimization. 

Figure  7-6a  gives  a montage  of  spatial  contours  of  dose 
distribution  on  the  coronal  plane  from  a modified 
cylindrical  dose  model  and  the  eight  optimum  arc  positions 
searched  by  the  Powell  method  for  case  2-2.  The  isodose 
shape  in  figure  7-8a  shows  a nearly  comparable  result  with 
that  obtained  from  the  exact  dose  model  near  the  crtical 
area  in  low  isodose  line  (Figure  7-8b) . The  dose 
distribution  with  the  final  search  values  in  Fig.  7-8b  gives 
a more  optimum  distribution  than  that  obtained  with  the 
starting  design  values  in  Fig.  7-Sc,  while  Keeping  the 
target  margin  and  minimizing  dose  near  the  critical  organs 
(lowered  from  25%  to  15%) . 


Isodosa  distribution  on  the  coronal  plane  from 
a modified  cylindrical  dose  model  (a)  and  exact 
dose  model  (b) , with  the  four  optimum  arc 
positions  searched.  The  isodose  lines 
displayed  are  from  10  to  90  in  increment  10  and 
normalized  by  maximum. 


Figure  7-B.  Isodose  distribution  on  the  coronal  plane  fron 
a modified  cylindrical  dose  model  (a) , and 
exact  dose  model  (b)  with  the  eight  optimum  arc 
positions  searched.  {c)  represents  the  isodose 
distribution  with  the  standard  arcs  and 
experimentally  determined  Isocenter  positions 
and  collimator  sizes.  The  isodose  lines 
displayed  are  from  10  to  SO  in  increment  10  and 
normalized  by  maximum. 


cylindrical  doe 


rotational  arcs,  which  are  generally  used  in  arc-based 
radiosurgery,  and  gives  dose  distributions  similar  to  those 
of  the  exact  dose  model  except  on  high  isodose  lines,  which 

to  the  exact  dose  model  and  much  more  efficient  in  finding 


critical  organ  through  computer-aided  design  optimization 
after  determining  optimum  isocenter  positions  and  collimator 
sizes  either  experimentally  or  by  computer-aided  design 
optimization. 


dose  distribution  to  the  target,  and  were  followed  by  an 
analytic  optimization  technique  to  search  for  the  best  arc 
parameters  to  further  minimize  the  dose  at  critical  points 


sizes  which  fit  the  target  shape  with  a high  isodose  surface 
surrounding  normal  structures.  A rule  based  search  method 


using  multiple  isocenters  and  multiple  collimators.  The 
second  step  is  to  find  optimum  arc  parameters,  with  which 
the  selected  critical  points  can  be  further  minimised  by 
varying  arc  parameters  while  maintaining  target  dose,  two 
previous  zero  order  search  algorithms  or  a new  rule-based 
search  method  could  be  used  with  the  cylindrical  dose  model. 
The  second  step  may  be  deleted,  if  further  minimization  to 
the  critical  points  is  not  necessary.  In  the  following 
examples,  we  consider  arbitrary  3-D  patient  data  which  can 
be  represented  by  a series  of  transvers  sectional  contours 
along  the  patient.  This  set  of  3-D  patient  data  contains 
target  cross-sections  and  relevant  normal  tissue  outlines  or 
critical  points. 

Problem  Chosen 

Case  3-1;  The  target  is  assumed  to  be  an  elongated 
cylindrical  shape  (height-2. 0 cm,  diameter-o.4  cm)  which  is 
represented  by  a series  of  contours  from  three  slices.  The 
critical  organs  are  also  located  at  two  different  positions, 
Cj{0,l,l)  and  Cj(0,-1,-1)  in  the  reference  BRW  coordinate 
system  (Fig.  7-9a) . 

case  3-2;  The  target  is  assumed  to  be  an  elongated 
cone  shape  (height-  2.4  cm,  diameter-  1.2  cm)  which  Is 
represented  by  a aeries  of  contours  from  four  slices  (Fig  7- 

We  optimize  case  3-1  and  case  3-2  using  an  analytical 
approach.  The  isocenter  positions  and  collimator  sizes  are 


critical  point 


dose  model  and  rule-based  search  methods  in  two  cases,  and 
is  fallowed  by  the  search  for  the  best  arc  parameters  (arc 
spacing  or  arc  weighting)  to  minimize  dose  at  the  critical 


Optimum  isocenter  positions  and  collimator  sizes 

The  isocenter  positions  and  collimator  sizes  are 

search  methods.  Our  aim  is  now  to  find  optimum  isocenter 

such  that  the  objective  function: 


is  maximized  subject  to 
Eqs.  7-2  and  3. 


snditions  described  in 


collimator  size 


standard  arcs) 


246 


Table  7-} QptinuM  Variabla  Vainaa 


Figure  7-10.  Isodose  distributions  on  the  plane  along  two 
isocenters  (a,c)  and  on  the  vertical  plane 
(b,d)  from  a exact  dose  model  with  the  optimum 
isocenters  and  collimator  sizes  searched  by 
the  rule-based  search  algorithms  for  case  3-1 
(a,b)  and  case  3-2  (c,d).  The  isodose  lines 
displayed  are  from  90  to  10  in  decrement  10 
and  normalized  by  maximum. 


approxinate  optimum  values  using  computer  search  methods  and 
approximate  dose  models,  it  is  useful  to  check  the  final 
dose  distribution  visually  using  the  exact  dose  model  with 
final  tuning  of  optimum  variables. 

The  isodose  shape  in  figures  7-lOa  and  d shows  a 
comparable  result  with  that  obtained  from  the  spherical  dose 
model  near  target  area  in  the  high  isodose  lines  (Figure 
7-lla  and  d) . The  dose  distribution  with  the  final  search 
values  in  Pig.  7-lOa  and  d give  a more  optimum  distribution 
than  that  obtained  with  the  starting  design  values  in  Pig. 
7-12a  and  d,  while  keeping  the  target  margin  and  maximizing 
dose  gradient.  It  is  known  that  the  dose  distribution  with 
a two  isocenter  approach  in  Pig.  7-lOa  and  d gives  a better 
distribution  than  that  obtained  with  a single  isocenter 
approach  in  Fig.  7-lia  and  d,  while  confining  the  target 
margin  within  the  70%  isodose  surface  saving  the  dose  to  the 
normal  structure  area. 

It  is  concluded  that  the  spherical  dose  model  is  useful 
in  calculating  dose  quickly  for  equally  spaced  arcs,  which 
are  generally  used  in  standard  stereotactic  radiosurgery, 
and  gives  dose  distributions  similar  to  those  of  the  exact 
dose  modal  except  at  low  isodose  lines,  which  makes  this 
approximate  dose  model  an  attractive  alternative  to  the 
exact  dose  model  and  much  more  efficient  in  finding  optimum 
isocenter  positions  and  collimator  sizes  to  fit  the  target 
shape  within  high  isodose  surface  while  maximizing  dose 


Isodose  distributions  on  the  plane  along  two 
isocenters  (a,c)  and  on  the  vertical  plane 
(b,d)  from  a spherical  dose  model  with  the 
optimum  Isocenters  and  collimator  sires 
searched  by  the  rule-based  search  algorithms 
for  case  3-1  (a,b)  and  case  3-2  (c,d).  The 
isodose  lines  displayed  are  from  90  to  10  in 
decrement  10  and  normalized  by  maximum. 


c 
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Figure  7-12.  Isodose  distributions  on  the  plane  along  two 
isocenters  {a,c)  and  on  the  vertical  plana 
(b,d)  form  a exact  dose  model  with  the 
starting  design  variables  (Isocenter  position 
and  collimator  siaa)  with  standard  arc 

isodose  lines  displayed  are  from  90  to  10  in 
decrement  10  and  normalized  by  maxinun. 


Figure  7-13.  Isodose  distributions  for  single  isooenter 

technique  with  a cylindrical  target  (a,b),  and 
cone-shaped  target  (c,d)  covered  by  80% 
isodose  surface.  The  isodose  lines  displayed 
are  from  90  to  10  in  decrement  10  and 
normalized  by  ntaximum . 


c 
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gradient  near  the  target  area  through  computer-aided  design 
optimization. 

Optimum  are  parameters 

The  isocenter  positions  and  collimator  sizes  were 
determined  to  shape  the  dose  distribution  to  the  target 
using  a spherical  dose  model  and  rule-based  search  methods 
in  case  3-1  and  case  3-2.  The  next  step  is  the  search  for 
the  best  arc  parameters  (arc  spacing  or  arc  weighting)  to 
further  minimize  dose  at  the  critical  points  chosen  in  case 

our  aim  is  now  to  find  optimum  arc  positions  (turntable 
orientation)  and  arc  weightings  for  the  individual  standard 
arcs  (100  deg  arc  length]  such  that  the  objective  function: 

S D,(»)  (7-la) 

is  minimized  subject  to  the  constraint  conditions  described 


d‘  < D,(x)  < 0?  Constraints  (7-2) 

x‘  S X S x“  side  constraints  (7-3) 


X = Design  vector  (turntable  orientations  and  arc 

weightings) 


D^(s)  = Doee  to  the  critical  organ  at  position  c 

Wj  =■  weighting  factor  to  evaluate  importance  of 
critical  organ  at  position  c 

d',  = Lower  and  upper  limits  of  the  desired  dose 

■ Lower  and  upper  limits  of  the  design 
variables 

The  c's  represent  positions  of  critical  organs  designated 
for  case  3-1.  is  relative  dose  to  the  critical  organ 
located  at  position  c,  and  is  represented  by  the  approximate 
dose  model  given  by  Eg.  5-91.  All  the  arc  weights  are 
assumed  to  be  one,  and  the  critical  organ  weights  are 
assumed  to  be  one.  The  criteria  of  uniformity  of  dose  Co 
the  target  region  is  also  defined  by  lower  and  upper  dose 
constraints  (e.g.  relative  dose  O.S  to  1.0)  at  preselected 
points  t's  within  the  target  area.  According  to  the 
physical  limitation  of  turntable  orientation  (e.g.  -90  to 
90),  the  lower  and  upper  limit  for  design  variables  are 
considered. 

Case  3-1  could  be  solved  using  eight  design  variables 
(i.e.  eight  turntable  orientations)  by  the  penalty  function 
approach  described  in  case  2-2.  The  another  simple 
alternative  might  be  to  utilize  a step  search  method  with 
constraint  conditions  described  in  Eqs.  7-2  and  3 even 
though  it  is  not  efficient  in  the  case  of  many  arc 
variables. 

Figures  7-14  give  montages  of  spatial  contours  of  dose 
distribution  on  the 


plane  along  from 


Figute  7-14.  Isodose  distribution  on  tha  coronal  plana 

aloncr  two  isocanters  from  a exact  dose  modai. 
The  optimum  isocenter  positions  and  coliimator 
sizes  were  searched  by  the  rule— based  search 
algorithm,  and  optimum  arc  positions  were 
searched  by  conjugated  direction  search 
algorithm. 


optinun  parameters  (optinun  Isocenters,  collimator  sizes  and 
algorithm  and  penalty  function  method  for  case  3-1.  The 
The  dose  distribution  with  the  final  search  values  in 


with  the  starting  design  values  in  Fig.  7-l2a  or  7-lOa, 
while  keeping  the  target  margin  and  minimizing  dose  near  the 
critical  points. 

It  is  concluded  that  approximate  dose  models  (spherical 
and  cylindrical)  are  useful  in  calculating  dose  quickly  for 

generally  used  in  standard  stereotactic  radiosurgery.  These 
approximate  dose  models  give  dose  distributions  similar  to 


design  optinizatior 
optimization  using 


ire  efficient  in  finding  optimum 
target  shape  within  a high  isodcse 


critical  organs  through  computer-aided 
Figure  7-15  shows  the  procedure  for 


aided  design  opciaizacion 


CHAPTER  8 
CONCLUSION 


distribution  for  arc-based  radiosurgery.  Tbe  experiaental 

technique  for  elongated  target  shapes.  Potential  studies 
for  shaping  2-D  (thin  plate)  or  3-D  (arbitrary)  targets  are 
expected  in  the  future.  The  slople  spherical  dose  nodal 

use  of  too  nany  isocenters  nay  not  be  desirable  to  shape  the 
conplicated  target  exactly,  since  It  gives  little  benefit 
with  nuch  increased  effort.  A conformal  therapy  using 
bean’s  eye  view  might  be  a better  approach  to  shape  the  more 
complicated  targets  with  an  computer-aided  design 
optimization  in  a future  study. 

we  conclude  that  in  optimal  dose  distribution  for 

forms  for  dose  calculation  in  3-D  space,  in  the  currant 
work  we  use  a spherical  or  cylindrical  dose  model  to 
simulate  the  exact  dose  distribution  for  a standard  arc,  and 
optimize  through  computer-aided  design  optimization.  Since 
the  dose  distribution  from  the  simple  dose  model  is  s 


jter-aided  design  optinizatlon  w 
1 efficient  and  practical  altern 


improvenents  are  expected  in  the  future.  One  is  finding  an 
arc  correction  factor  for  any  length  of  arc  in  addition  to 


It  possible  to  use  the  start  and  stop  positions  of  gantry 
angle  as  variables  in  CAD  optinization.  Since  the  dose 


head  contour  vary,  it  nay  not  be  necessary  to  correct  tor 
different  target  positions  and  head  contours.  In  the 

It  must  be  interesting  points  to  consider  nultlple  targets 

situations.  Another  tash  is  finding  better  approximate 
forms  using  a different  analytical  form. 

We  may  consider  the  possibility  of  using  computer- 
using rule-based  or  step  search  method.  However,  the  use  of 


efficient  optimization  criteria  and  that  of  specifying  the 
formalisms,  it  may  be  desirable  to  develop  a feu  standard 


optimal  result  by  varying  some  parameters  using  multiple 
steps.  The  CAD  optimization  technique  should  be  initiated 

treatment  planner  and  the  three*dimensional  treatment 
planning  system.  Another  task  is  to  develop  a method  which 
simulates  a real  patient  treatment  with  a reference  patient 


oriteria.  The  statistical  approach  to  optimization 

probabilities,  and  normal  tissue  complication  probabilities 

all  major  factors  or  statistical  information  can  be 

optimization  criteria  might  be  good  to  implement  CAD 
optimization  techniques,  and  could  be  used  until  all  the 
biologic  information  are  obtained.  The  selection  criteria 

fxuictlon  and  constraints  are  very  important  to  the  success 
of  the  computer  analytic  optimization.  The  simplified  dose 

analytical  optimization  techniques. 

"Hooke  and  Jeeves  method"  and  "Powell  method"  were  used  with 


an  interior  penalty  function  approach.  Although  these 
methods  showed  some  good  results  in  simple  examples,  this 
method  may  not  be  the  best  choice.  After  testing  some  other 
optimization  methods  (zero  order,  first  order,  or  seoond 
order  with  penalty  approach  or  direct  search  method),  the 
most  reliable  method  nay  be  applied  to  practical 
optimization  problems. 
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